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SUMMARY 


An extensive literature review on the effect of interface/ 
interphase on composite properties was followed by experiments 
which illustrate the significance of interphase in composite 
performance. 

The review covered adhesion theory, wetting characteristics 
of carbon fiber, load transfer mechanisms in a fibrous composite, 
test methods to measure and evaluate fiber-matrix interfacial 
bond strengths, effect of moisture at the interface and properties 
of the interface/interphase, including impact toughness, on com- 
posite performance. A critical discussion of the literature is 
presented along with appropriate citations and references. 

In contrast to adhesion in glass fiber /epoxy composites, 
carbon fiber/epoxy adhesion is not simply promoted by "coupling 
agents" but rather it results from physical and/or chemical inter- 
actions with the matrix. Much work has been done to characterize 
the surface energetics of carbon fiber but no correlation has 
been identified between wetting behavior and interfacial bond 
strength. Many techniques have been developed to study fiber/ 
matrix interfacial bond strength but few are applicable to car- 
bon fiber. Preliminary data have shown that single filament 
critical length determination is applicable to evaluate carbon 
fiber - epoxy bond strength. This technique was adapted in this 
study to evaluate finish variants on carbon fiber in epoxy com- 
posites . 

Moisture at the interface of carbon fiber/epoxy composites 
has been shown to have little effect on composite properties. 

The decrease in composite properties in humid environments is 
primarily due to plasticization of the matrix by water. 

Improvement of impact toughness of carbon fiber/epoxy com- 
posites can be achieved by modification at the fiber/matrix inter- 
face. However, this is usually associated with a significant 
decrease in mechanical properties. This is because weak inter- 
facial bonds allow more impact energy to be absorbed, thereby 
increasing the impact toughness. A few workers have suggested 
that "tailoring" interphase properties can improve toughness 
without lowering the mechanical properties, but, this concept 
has not yet been demonstrated experimentally. 

In this study, experimental results were obtained which show 
the effect of interphase on composite performance and illustrate 
its significance. Various finish variants were formulated, based 
on different chemical and mechanical properties, and applied to 
Celion® 6000 carbon fiber. Wetting behavior, studied by the 
Wilhelmy wetting force technique, and interfacial bond strength, 
evaluated by single filament critical length determination were 
characterized for these fibers. Composite properties at room and 
elevated temperatures and impact toughness were also measured. 
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It is concluded that good wetting is essential to provide 
good adhesion but better wetting properties do not necessarily 
translate to better bond strength. Also, it is shown that with 
a weak interfacial bond strength, impact toughness is improved. 
But, most significantly, a 30% increase in impact toughness was 
measured without loss of room temperature and elevated tempera- 
ture composite mechanical properties by "tailoring" the inter- 
phase properties. 
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INTRODUCTION 


In the last decade, much attention has been focused on fiber- 
reinforced polymeric matrix composites. These composite materials 
are both versatile and complex. Their versatility stems from a 
wide choice of constituent materials available and from the variety 
of ways in which composites can be fabricated to provide a combina- 
tion of desired properties. These tailor-made properties are fre- 
quently not available and cannot be achieved in conventional isotro- 
pic materials. On the other hand, they are complex by virtue of 
their chemical and mechanical nature and thus should not be regarded 
as a single material, but as a material system. Inherent in this 
multi-phase material system is the fiber/matrix interface which is 
the critical link that provides the structural integrity for the 
fibrous composites. 

Conventionally, the interface is regarded as the bond between 
fiber and matrix. This interfacial bond is considered to have zero 
thickness and to result from the interaction between the fiber surface 
and the matrix material. The interphase region of a fibrous composite 
is the area immediately adjacent to and including the interface and 
extends a finite distance into the bulk matrix material. This concept 
is illustrated in Figure 1. 

An extensive review of the literature indicates that much 
work has been done on fiber/matrix interface characterization. 

These include the evaluation of fiber surface energetics via wet- 
ting force measurement; surface elemental analysis via ESCA, 

AUGER, XPS, etc. and surface topography /morphology studies via 
SEM. However, little effort has been spent to correlate this 
information with composite performance. In short, the role of 
the interface in determining composite properties is not known. 
Furthermore, little attention has been given to the understand- 
ing of the nature of the interphase and how it affects composite 
properties. In particular, the load transfer mechanism in the 
interphase region is not clear. This report is intended to 
clarify these areas of uncertainty and to illustrate the influ- 
ence of the interphase region on composite performance. 
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BACKGROUND AND LITERATURE REVIEW 

ROLE OF THE INTERFACE IN COMPOSITE MATERIALS 

Among many factors that govern the characteristics of a 
fibrous composite, it is clear that the adhesion at the fiber/ 
matrix interface plays an important role in structural integrity 
of these two-phase materials . Since there are a limited number 
of ways to directly study adhesion at the interface, composite 
failure mode and shear strength are commonly employed to assess 
interfacial bond strengths. 

Two extreme cases generally occur in interfacial failure 
mode studies. For weak interfacial bonding, fiber "pull out" 
is responsible for the rupture. This is because the fibers are 
not bonded by the matrix and they slide in their cavities. 
Therefore, load cannot be transferred in such a discontinuous 
medium. In this case, the composite fails in a shear mode and the 
shear strength is usually low. On the other hand, if the inter- 
facial bonding is strong, then the matrix will transfer the load 
to the fibers until it reaches the breaking strength of the 
fibers. In this case, the failure is sudden and catastrophic. 

A higher shear strength is associated with this kind of tensile 
failure. 

Virgin carbon fibers exhibit lower bondability with polymer 
matrixes (e.g., epoxy, polyester) and surface treatments are 
employed to increase the interfacial bond strength. All commer- 
cially available carbon fibers are surface treated. These surface 
treatments are generally oxidative in nature and provide 
functional groups on the fiber surface which increase wettability 
and/or bondability. Much of the early work has emphasized 
increases of composite shear strength. 

Goan et al(!) studied a sodium chlorate/sulfuric acid 
oxidation process which doubled composite short beam shear 
strength. However, this process degraded the fiber tensile 
strengths. A 4% loss in fiber strength (impregnated strand 
tensile test) was observed. Druin et al< 2 ) studied gas phase 
oxidation processes where similar results were obtained. 

Carbon fiber strength is surface flaw sensitive. Extensive 
oxidation/etching can increase surface flaws which leads to a 
reduction in properties. Ehrburger et al, Molleyre and 
Bastick, ( 4 ) Fitzer et al( 5 ) and Duffy( 6 ) studied various types 
of oxidation processes to optimize the incorporation of functional 
groups onto the fiber surface and reduce the possibility of exten- 
sive oxidation/etching. 

Ehrburger et al( 7 » 8 ) found that, by the insertion of weak 
acidic groups, interlaminar shear strength (ILSS) could be en- 
hanced just as much as with the incorporation of strong acidic 
groups on AC (type II) , ex— acrylic carbon fiber surfaces. However, 
the mode of failure was different. A shearing mode of fracture was 
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observed with weak acidic groups while tensile failure was ob- 
served with strong acidic groups. A possible explanation suggested 
by the author for this observation was that the weak acidic groups 
increased the surface energetics of the fibers and resulted in 
better wetting and good adsorption which enhanced the coupling. 

The strong acidic groups - mainly carboxylic or phenolic groups - 
are coupled with the matrix via tighter and stronger primary 
bonds at the interface which create a brittle composite. 

Although little is known about the adhesion mechanisms 
involved in a fibrous composite material, it is clear that the 
role of the interface in a composite material is significant. 
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EFFECT OF INTERPHASE PROPERTIES ON COMPOSITE PERFORMANCE 

Commercially available carbon fibers are usually sized with 
an epoxy compatible finish to enhance handling characteristics. 
Epoxy-based systems are generally used since epoxy resins are com- 
monly used as components of the matrix materials and such mate- 
rials can be easily applied to the fibers. These sizes also act 
as lubricants as well as "loose" binders which prevent fiber dam- 
age and minimize fiber fuzz. The loose binder should not hinder 
the spreading of the fibers during the prepregging process, but 
should allow the resin to penetrate bundles and wet individual 
fibers to eliminate dry spots or voids. 

Sizing can be accomplished by passing tows of fibers through 
a bath containing a dilute solution of the finish material. Fol- 
lowing such a coating, the fibers must be dried to remove solvent. 
Better adhesion properties were observed by Drzal et al(l) and 
Goan et al(2)with commercial epoxy-sized fibers. This phenomenon 
can be explained on the basis that the fibers were sized with a 
low viscosity fluid, thus enhancing wetting. Furthermore, the 
spontaneous spreading of the sizing solution minimizes trapped 
air pockets at the interface between fiber and finish. 

The overall effect of size materials on composite properties 
is not fully understood. This is partly because the size materials 
are applied onto the fiber in very thin layers, normally 1.5% by 
weight of the fiber, which is equivalent to an approximately 800 A 
thick coating (Appendix A) . 

Few analytical techniques can be applied to study the inter- 
phase of a carbon fiber epoxy composite because of its small 
dimension. Furthermore, the fiber geometry, and the low atomic 
number elements on the surface further limit the applicable analyti- 
cal techniques. These add to the difficulties in understanding the 
role of the interphase in composite properties. However, some of 
the effects of interphase can be observed in mechanical tests. 

Arai et al(3) showed an increase in both interlaminar shear 
strength (ILSS) and flexural strength of a carbon fiber epoxy com- 
posite by electrolytic polymerization of vinyl monomers onto the 
carbon fibers. Pinchin and Woodhams(^) pyrolyzed graphite onto car- 
bon fiber surfaces and increased the surface energetics. This was 
reflected in an increase in shear (ILSS) and flexural strengths. 
Furthermore, a slight improvement on impact toughness was also 
measured. An hypothesis, based upon this phenomenon, is that the 
pyrolytic graphite deposit forms a sheath around the fiber which 
creates an additional interface. During the debonding and/or crack 
propagating processes, this additional interface absorbs more energy 
by exposing additional free surfaces. 

Broutman and Agarwal^ • 7 ) did a theoretical study on composite 
properties with an interphase layer. Assuming both phases behave 
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elastically with perfect bonding at the interface, it was shown 
that composite strength reached a maximum when the interphase 
modulus reached a value of 69 MPa (10 4 psi) . It should be noted 
that Broutman used the term interface modulus but, according to 
the definitions of interface and interphase described earlier, 
interphase modulus is more appropriate because interface is re- 
garded as zero thickness . Further increases in the interphase 
modulus do not change the composite strength. This is shown in 
Figure 2. However, composite modulus continues to increase with 
increasing interphase modulus (Figure 3) . This further increase 
in modulus without composite strength changing further, indicates 
that composite elongation can be maximized at this critical inter- 
phase modulus (Figure 4) . 

Although Broutman' s theoretical study was not supported by 
experimental data, the critical interphase modulus determined in 
his study was significantly lower than the modulus of the matrix. 

This is similar to the flexible interlayer concept discussed by 
Benedetto and Nicolais (5). This concept attributes the effective- 
ness of finishes to stress relief in the interphase region. Since 
there is a difference in thermal expansion coefficient of the fiber 
and the matrix, significant residual stress will accumulate at the 
interface due to elevated temperature curing as well as shrinkage 
of the matrix during polymerization. However, by incorporation of 
a flexible interlayer, residual stresses can be relieved. Although 
these two concepts are different in approach, they both conclude 
that the properties of the composite can be improved through selec- 
tive tailoring of the interphase. 

. (12) 

Marom and Arridge studied the stress patterns of soft inter- 

layers on stainless steel inclusions. Reduction of stress concen- 
tration at the interface was observed and an improvement in the ul- 
timate tensile strength was measured. However, Marom stated that 
optimization of this soft interlayer is necessary to prevent poor 
transverse properties. This effect was reported by Gatward and 
Hull, (13 ) when they coated a flexibilized polyvinyl butyral polymer 
onto glass fiber and a decrease in transverse tensile properties 
was measured. 

Taniguchi et al ( 9 ) incorporated a flexibilized epoxy at the in- 
terphase of a carbon fiber/epoxy composite and preliminary data 
showed improvement in flexural and shear strength. Kardoes et al(10>H) 
and Lavengood et al(14) coated boron fibers and glass fibers with a 
low cross-linking density epoxy system to a thick, flexible inner- 
layer. Increase in mechanical properties, both moisture conditioned 
and unconditioned, were measured. A 1000% increase in torsional 
fatigue life was measured and Figure 5 shows the improvement of 
transverse properties of glass fiber/epoxy composites by incorpora- 
tion of a soft interphase. 
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COMPOSITE STRENGTH AS A FUNCTION OF 
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COMPOSITE MODULUS AS A FUNCTION OF 
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*From Broutman and Agarwal (Ref. 7) 
Reproduced by permission. 
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FIGURE 5. TRANSVERSE STRENGTH OF GLASS FIBER REINFORCED 
COMPOSITES BEFORE AND AFTER TWO HOURS IN 
BOILING WATER* 


*From Lavengood and Micheno (Ref. 14) . 
Reproduced by permission. 
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( 8 ) 

Riess et al ' developed an interphase material with poly- 
Cisoprene-b- (styrene-alt-maleic anhydride)] block copolymer: 



(Isoprene) (SAM) 


According to the compatibility of the polymers, the polyiso- 
prene will agglomerate to form micellar structures at the inter- 
phase while the SAM block is mixed with the matrix resin to form 
a continuous phase. This inclusion of a second, elastomeric phase/ 
is to relieve residual stress and to damp the elastic energy as 
discussed earlier. A proposed interphase structure is shown in 
Figure 6 . 




FIGURE 6. MICELLE STRUCTURE OF THE COPOLYMER 
MICELLAR ADSORPTION* 

*From Riess et al (Ref. 8). 

Reproduced by permission 


This PI -SAM interphase material can couple to the carbon 
fiber surface through the dipole-dipole interaction sketched in 
Figure 7 . 
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§ elastomer 

SAM 

MATRIX OF COMPOSITE 


FIGURE 7. FIXATION OF BLOCK COPOLYMERS ON CARBON 
FIBRES BY DIPOLE-DIPOLE BONDING* 


The introduction of this "soft" interphase not only increased 
the shear strength of the composite, both unconditioned and con- 
ditioned in water at 50°C, but also improved the impact toughness. 
Some of these results are shown in Figures 8 and 9 . 



FIGURE 8. SHEAR STRENGTH OF COMPOSITES AFTER WATER TREATMENT 
AT 50 °C* 

© RESIN 
© RESIN + FIBRES 
© RESIN +GRAFTED FIBRES 

*From Riess et al. (Ref. 8) . 

Reproduced by permission. 


14 




FIGURE 9. SHEAR STRENGTH AND RESILIENCE OF: 
RESIN (EPOXY) 

RESIN + UNTREATED FIBRES 
RESIN + GRAFTED FIBRES* 


*From Riess et al (Ref. 8) . 
Reproduced by permission. 


In reviewing the theoretical concepts together with the 
limited experimental data, a soft interphase between fiber and 
matrix clearly improves some of the composite properties. The 
objective of this report is to provide further understanding of 
the role of interphase in composite properties and to define 
guidelines and principles for the selection of optimal interphase 
materials to improve overall performance. 
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ROLE OF INTERFACE/INTERPHASE IN COMPOSITE TOUGHNESS 

Toughness in composite materials is not a well defined 
property. This is due in part to the difficulties in designing 
test methods to study the toughness of a multi-phased aniso- 
tropic material. In general, toughness is measured by the 
resistance of the material to crack growth and propagation. 
However, for multi-phase anisotropic materials, this is not 
totally representative. 

An increasingly common method in toughness measurement is 
instrumented impact. This method utilizes a load cell mounted 
onto an impacter. The load cell is capable of recording load 
during impact in microseconds. A microprocessor then outputs a 
load deflection curve together with the energy absorption 
profile during the course of impact. The maximum force on the 
load deflection curve is taken as the required load for penetra- 
tion . 


Due to the multi-phase characteristic of a composite, its 
toughness can be improved in several ways. A toughened matrix 
and/or a higher strain fiber are both capable of increasing the 
toughness of a composite. Although little is known about the 
interface and interphase, it has the potential to significantly 
influence the toughness of composites. 

The interface can affect the toughness of a composite by 
providing various energy absorbing mechanisms like debonding, 
fiber stress relaxation and fiber pullout during fracture. 

In the case of weak interfacial bonds, debonding takes place in 
advance of the crack tip. The degree of debonding, as well as 
the energy absorbed due to the debonding process, depends largely 
upon the inter facial bond strength. In the opposite case of 
strong interfacial bonds, fibers will break prior to debonding 
and energy is absorbed due to stress relaxation at the ends of the 
broken fibers. Furthermore, energy is also absorbed due to some 
fiber pullout which is related to the strength of interfacial 
bonds. When all of these energies are accounted for, it is 
expected that the weak interfacial bond will provide very high 
fracture toughness in composites. This phenomenon has been 
observed by many workers. (1“4) 

However, by tailoring the interphase properties, a signifi- 
cant improvement in composite toughness can be obtained without 
a severe loss of mechanical properties due to weak interfacial 
bond strength. Broutman et al(5-7) have analytically demonstrated 
that interphase properties can be tailored to obtain both toughness 
and strength. In general, a softer material at the interphase is 
preferred. 

Plueddemann (8) incorporated a rubber interphase in glass/ 
epoxy composites. Improved toughness together with relatively 
good adhesion strength was obtained. Hancox and Wells (9) coated 
different amounts of silicone rubber at the surface of carbon 
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fibers. A 107% increase in work of fracture without measurable 
decrease in flexural strength was determined with 1% by weight 
silicone rubber at the interphase of a carbon fiber/epoxy com- 
posite. Williams (10) polymerized polyvinylalcohol and polysul- 
fone onto carbon fiber and enhancement in shear strength, notched 
fracture toughness and low-velocity impact toughness were achieved. 
Subramanian and Jakubowski (11) coated carbon fibers with various 
kinds of polymers by electropolymerization and demonstrated that 
impact toughness and shear strength can simultaneously increase 
when optimum interphase properties were incorporated. 

In addition to the mechanical properties of the interphase, 
the physical thickness of the interphase is another critical 
parameter that needs to be optimized in order to achieve the 
ultimate properties of composite materials. As discussed by Han- 
cox and Wells, (9) when silicone rubber exceeded the 1% by weight 
level, a dramatic decrease in flexural strength was observed. 

Similar observations were also reported by Cavano and Winters (12) 
and Peiffer. (13) 

Another approach was taken by Atkins (I 4 ) where intermittent 
bonding along the fiber surface was suggested. The weakly bonded 
sections allow debonding and fiber pullout during fracture while 
the strongly bonded areas enable load transfer. The combination 
of these could then simultaneously provide high fracture toughness 
and good tensile strength. Furthermore, Reiss et al(15) incorpor- 
ated a two-phase block copolymer, polyisoprene-polystyrene-maleic 
anhydride, at the interphase. Where the polyisoprene is providing 
the energy absorbing characteristic and the stryene-maleic anhy- 
dride block (SAM) is enhancing the compatibility and joint strength 
of the fiber with the matrix, increases in shear strength and 
toughness were both measured. In addition, this tailored interphase 
enhanced the moisture-conditioned shear strength. 

In view of the results discussed above, composite toughness 
can be increased without accompanying degradation of the remain- 
ing mechanical properties. This can be achieved by proper selec- 
tion of interphase properties, such as transition temperature, 
modulus, strength, etc. and the physical characteristics such as 
the quantity and thickness of this interphase. 
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EFFECT OF MOISTURE ON THE CARBON FIBER/EPOXY INTERFACE AND INTERPHASE 

Moisture is usually found on all solid surfaces and water 
molecules are usually the inhibitors/contaminants for adhesive 
bonding. Water is generally attracted end adsorbed onto the sur- 
face of solids through hydrogen bonding to the surface hydroxyl 
groups . The water molecules are normally one or two monolayers 
deep on the solid surface but this is enough to induce stress 
corrosion and assist crack propagation. 

Evidence of adsorbed water molecules at glass fiber/epoxy 
interfaces has been shown and its detrimental effect on composite 
properties were discussed by Bascoro. (1) According to Bascom, 
water molecules are attracted by the highly concentrated surface 
hydroxyl groups on the glass fibers and the reduction in composite 
properties are due to stress corrosion at the interface as well 
as in the bulk matrix. 

Although carbon fibers also contain many functional groups 
on their surfaces, moisture has a smaller effect on carbon/epoxy 
composites. According to Gauchel and Nash, (2) only minor reduc- 
tions in composite properties are observed after long-term exposure 
to water (12 months) . Furthermore, these property reductions 
were primarily due to the plasticization effect of water on the 
epoxy matrix. Springer and Loos, (3) and Augl(4,5) showed that 
by characterizing effects of moisture, on epoxy matrix properties, 
the moisture effect on carbon/epoxy composites’ performance can be 
predicted with a high degree of confidence. In other words, com- 
posite performance reduction is primarily due to bulk matrix 
plasticization . 

In addition, McKague et al(6) at General Dynamics have studied 
the moisture effect on carbon/epoxy composites as a function of 
the sizing materials . McKague concluded that the size variants 
studied showed no effect on total moisture absorptivity (Figures 10 
through 15) or glass transition temperature (Table I) , nor did it 
affect the resistance to thermal cycling between — 54°C and 124°C 
( -65 °F and 255°F), (Table II). 

In conclusion, for carbon fiber/epoxy composites, the decrease 
in shear strength after moisture conditioning is primarily due to 
plasticization of the matrix by water. This is contrary to glass 
fiber/epoxy composites, wherein irreversible damage has been re- 
ported due to moisture degradation of chemical bonds at the inter- 
face. This has not been shown for carbon fiber/epoxy composites. 
Although moisture may have some effect on the interphase region in 
carbon fiber composites, experimental techniques are not available 
which can separate the effects of plasticization of the interphase 
from plasticization of the bulk of the matrix itself. Furthermore, 
there is no evidence of chemical degradation of the interphase region 
by moisture at ordinary use temperatures. Therefore, no experimental 
effort will be devoted in this study to demonstrate the significance 
of moisture at the interface and interphase. 
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FIGURE 10. MOISTURE ABSORPTION OF CELANESE FURNISHED CELION®6000/5200 

LAMINATES WITH AND WITHOUT NOMINAL SIZING FINISH* 
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*From McKague et al (Ref. 6) 

Reproduced by permission. 
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FIGURE 11. MOISTURE ABSORPTION OF CELANESE FURNISHED CELION®6000/5200 
LAMINATES— SERIES 743-A AND 115 -A SIZING* 


*From McKague et al (Ref. 6) 
Reproduced by permission. 
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FIGURE 12 . MOISTURE ABSORPTION OF CELANESE FURNISHED C ELION®6000/5200 
LAMINATES —SERIES 743-A AND 101-A SIZING* 
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*From McKague et al (Ref. 6) 
Reproduced by permission. 
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FIGURE 13. MOISTURE ABSORPTION OF CELANESE FURNISHED CELION®6 000/5200 
LAMINATES --SERIES 743-A AND 747-A* 


*From McKague et al (Ref. 6) 
Reproduced by permission. 
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FIGURE 14. MOISTURE ABSORPTION OF CELANESE FURNISHED CELION®6000/5200 
LAMINATES — SERIES 743-A AND 120-A* 
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*From McKague et al (Ref. 6) 
Reproduced by permission. 




<20 NIHS 1H3I3H 



EXPOSURE TIME (50. RT. OF DRY5) 


FIGURE 15. MOISTURE ABSORPTION OF CELANESE FURNISHED CELION®6000/5200 
LAMINATES— SERIES 743-A, 745-A AND 747-A* 


*From McKague et al (Ref. 6) 
Reproduced by permission. 
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TABLE I. GLASS TRANSITION TEMPERATURES 
ON INTERPHASE VARIANTS* 


Sample 

Dry Tg °C 

(°F) 

Moisture (%) 

Wet Tq °C 

(°F) 

No Finish 

232 

(450) 

0.70 

166 

(330) 

Normal 1.1% 
Finish 

232 

(450) 

0.67 

166 

(330) 

0.5% 5210 
Epoxy 

232 

(450) 

0.80 

160 

(320) 

1% 5210 
Epoxy 

232 

(450) 

0.76 

168 

(335) 

2% 5210 
Epoxy 

232 

(450) 

0 . 78 

166 

(330) 

1% 5210 
Epoxy 

(emulsion) 

232 

(450) 

0.73 

160 

(320) 


Note: Each value is the average of two specimens. 

Maximum pair difference was 43.3°C (110°F) . 


* From McKague et al (Ref. 6) . 
Reproduced by permission. 
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TABLE II 


THERMAL CYCLING DATA* 


0° Flexural Properties 
Strength Modulus 



Uncycled 

Cycled 

Uncycled 

Cycled 

Sample 

MPa 

ksi 

MPa 

ksi 

GPa 

Msi 

GPa 

Msi 

No Finish 

1931 

280 

1924 

279 

142.1 

20.6 

144.8 

21.0 

Normal 1.1% 
Finish 

1869 

271 

1952 

283 

140.7 

20.4 

144.8 

21.0 

0.5% 5210 
Epoxy 

1876 

272 

1897 

275 

138.6 

20.1 

136.6 

19.8 

1% 5210 
Epoxy 

1897 

275 

1938 

281 

135.2 

19.6 

138.6 

20.1 

2% 5210 
Epoxy 

1841 

267 

1814 

263 

131.0 

19.0 

132.4 

19.2 

1% 5210 
Epoxy 

1938 

281 

1986 

288 

142.1 

20.6 

139.3 

20.2 


(Emulsion) 


*From McKague et al. (Ref. 6). 
Reproduced by permission. 
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ADHESION THEORY 


Adhesion, in physical chemistry, is the force of attraction 
between a solid surface and a second phase. The second phase can 
be a solid or a liquid. Adhesion leads to sorption, which may be 
adsorption on a surface or adsorption into a surface layer. Ad- 
hesion can be due to electrostatic forces, to van der Waals forces 
or to "chemical" valence forces. In the latter case, one speaks of 
"chemisorption." The magnitude of the adhesive forces is derived 
from equilibrium measurements of adsorption energies by thermodynamic 
methods . 

In adhesion technology, only the interaction between a solid 
surface and a second liquid or solid phase is termed adhesion. The 
technical process of producing adhesion between two solids is called 
adhesive bonding. The bonded system is not necessary in a thermo- 
dynamic equilibrium; frequently the bonding process is irreversible. 
Therefore, surface interactions cannot be separated experimentally 
from bulk properties of the system. 

In addition, the state of fundamental knowledge concerning the 
nature of the forces of attraction between bodies, determination of 
the magnitude of such intrinsic forces, and their relation to measured 
adhesive joint strengths undoubtedly lags behind state-of-the-art 
applied technology. One of the main reasons why developing theory 
has followed behind technology is that the science of adhesion is a 
truly multi-disciplined subject. The adhesion scientist often needs 
to consider aspects of surface chemistry and physics, chemistry, 
polymer chemistry, polymer physics, rheology, stress analysis and 
fracture phenomena to interpret his data fully. 

The mechanisms of adhesion are still not fully understood and 
many theories are to be found in the current literature. Much of 
present confusion results from the fact that test methods commonly 
employed to measure the strengths of adhesive joints are not well 
suited to theoretical analysis. They introduce geometrical factors 
and loading factors which are difficult to analyze, and the measured 
joint strength includes indeterminate contributions from rheological 
energy losses in the adhesives and substrate. Thus, although the 
intrinsic adhesion forces acting across the adhesive/substrate inter- 
face may affect joint strength they are usually completely obscured 
by other contributions, and information concerning the magnitude 
of such forces may only be indirectly obtained. This inability to 
measure the interfacial interactions has been the main obstacle to 
the development of a comprehensive theory of adhesion. 
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The four principal mechanisms of adhesion which have been pro- 
posed and widely discussed are: 

1. Mechanical Interlocking 

2. Diffusion Theory 

3. Electronic Theory 

4 . Adsorption Theory 

These are discussed in the following paragraphs. 

Mechanical Interlocking 

This theory proposes that mechanical keying or interlocking of 
the adhesive into the irregularities of the substrate surface is 
the major source of intrinsic adhesion. Much work^i - ^) i n this area 
has shown the importance of surface topography/rugosity in intrinsic 
adhesion forces. However, the observed increase in measured joint 
strength with increasing surface rugosity may be attributable to 
other mechanisms. 

Diffusion Theory 


The diffusion theory of adhesion is based upon the intrinsic 
adhesion of high polymers to themselves (autohesion) , and to each 
other, due to the mutual diffusion of polymer molecules across the 
interface. This requires that the macromolecules or chain segments 
of the polymers (adhesive and substrate) possess sufficient mobility 
and are mutually soluble, that is, they possess similar values of 
the solubility parameter. The solubility parameter, <5, is defined 

by 


5 


AH v - RT 


V 


( 1 ) 


where 

AH v = molar heat of vaporization 
R = gas constant 

T = temperature (°K) 

V = molar volume 

Voyutskii(8) and others (9-13) have presented the evidence of 

molecular diffusion and their significant contributions to the in- 
trinsic adhesion. However, with polymer systems which are highly 
crosslinked and have a high degree of crystallinity, diffusion 
theory is not likely to be the principal mechanism of adhesion. 
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Electronic Theory 


The electronic theory of adhesion arises from the different 
electronic bond structures of the adhesive and the substrate which 
cause some electron transfer on contact to balance Fermi levels. 

This will result in the formation of a double layer of electrical 
charge at the interface which contributes to the intrinsic adhesion. 
This theory was proposed by Deryaguin d 4 ”16 ) an $ supported by other 
workers. ( 17 24 ' However, the arguments were not strong enough to 
show its significant contribution to the intrinsic adhesion. None- 
theless, surface electronic states should not be overlooked. 

Adsorption Theory 

The adsorption theory of adhesion is by far the most generally 
accepted theory and has been discussed in depth by many authors.( 2 ^“ 28 ) _ 
This theory proposes, provided that sufficiently intimate inter- 
molecular contact is achieved at the interface, that materials will 
adhere because of the surface forces acting between the atoms in the 
two surfaces . The most common of such forces are the van der Waals 
forces (secondary bonds) . In addition, chemisorption may well occur 
and thus ionic, covalent and metallic bonds may operate across the 
interface (primary bonds) . 

The thermodynamic work of adhesion required to separate a unit 
area of two phases forming an interface, Wa, may be related to the 
surface free energies by the Young-Dupre equation. ( 2 ^) i n the absence 
of chemisorption and inter-diffusion, the reversible work of adhesion, 
Wa/ in an inert medium may be expressed by 

W A =Y 1 +Y 2“ Y 12 ( 2 ) 

where 

Y 1 = surface free energy of phase 1 
Y 2 = surface free energy of phase 2 
Y 12 .= interfacial free energy between the phases 

This equation only applies to a solid/liquid or liquid/liquid inter- 
face. However, by assuming that the surface energy of a liquid 
does not change significantly where it solidifies isothermally and 
ignoring the shrinking stresses, it may be applied to solid adhesives/ 
substrate interfaces. 

Dukes (30) an £ others (31-34) have been trying to correlate 
measured strengths of adhesive joints with the thermodynamic work 
of adhesion, Wa* However, only some of the authors have confirmed 
correlations while others have found no correlations at all. The 
major cause of this conflict arises from the use of test methods for 
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measuring the strengths of adhesive joints which are not well 
suited to theoretical analysis. Moreover, correlation can only 
be expected when the locus of joint failure is purely inter facial. 
Therefore, even though the thermodynamic work of adhesion affects 
joint strength, it is often obscured by other contributions. 

Although it has been shown that secondary bonding forces alone 
can give rise to acceptable joint strengths at the interface, many 
adhesion scientists believe that the additional presence of primary 
bonding can increase the joint strength and is certainly a must 
for establishing environmentally stable interfaces. Evidence of 
the presence of primary interfacialy bondinq has been shown via 
the use of sophisticated, surface-specific, analytical tech- 
niques such as Laser-Raman spectroscopy, (35) x-ray photoelectron 
spectroscopy , ^ 36 ) secondary-ion mass spectroscopy , (37 , 38) an( j 
attenuated total reflectance infra-red spectroscopy. ( 39 ) 

It was recognized by adhesion scientists that the establish- 
ment of an intimate molecular contact at the interface is a necessary, 
though sometimes insufficient, requirement for developing strong ad- 
hesive joints. In other words, the adhesive needs to be able to be 
spread over the solid substrate surface and needs to displace air 
and any other contaminants that may be present on the surface. 

Ideally, an adhesive should conform to the following conditions: (40) 

1. Adhesives in the liquid state must exhibit a zero or 
near-zero contact angle with the solid substrate. 

2. The adhesive should have a relatively low viscosity (a 
few centipoises) at some time during the bonding opera- 
tion. 

3 . The adhesive should be brought together with the sub- 
strate at a rate and in a manner that will promote the 
displacement of any trapped air at the interface. 

In order to assess the ability of a given adhesive/substrate 
combination to meet these criteria, it is necessary to consider 
wetting equilibria, to ascertain values of the surface free energies 
of the adhesive and substrate and the free energy of the adhesive/ 
substrate interface. In addition, the kinetics of the wetting 
process should also be examined. 
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WETTING PHENOMENA 


Surface tension is a direct measurement of intermolecular 
forces . The tension in surface layers is due to the attraction 
forces exerted on the surface molecules by the bulk material. This 
internal attraction force leads to a reduced number of molecules at 
the surface which result in an increased intermolecular distance . 

This increase in intermolecular distance requires external work to 
be done which explains why surface tension/surface free energy exist. 

The most common type of physical surface attractive forces are 
van der Waals forces which can be attributed to: (a) dispersion 

(London-d) forces arising from internal electron motions which are 
independent of dipole moments and; (b) polar (Keesom-p) forces aris- 
ing from the orientation of permanent electric dipoles and the in- 
duction effect of permanent dipoles or polarizable molecules. The 
dispersion forces are usually weaker than the polar forces but they 
are universal and all materials exhibit them. 

Another type of force that may operate is the hydrogen bond. 
Hydrogen bonding is the attraction between a hydrogen atom and a 
second, small and strongly electronegative atom such as fluorine, 
oxygen or nitrogen. A secondary bond is actually formed between 
the hydrogen atom and the electronegative (host) atom through an 
electron-pair on the host species. 

Equilibrium wetting occurs when the liquid and solid phases 
are in contact and can no longer be perturbed by external forces. 

The equilibrium condition for the wetting of a solid surface by a 
liquid is a three-phase equilibrium of solid, liquid and vapor. As 
shown in Figure 16, one can describe the equilibrium point of 


Vapor y 

LV 



FIGURE 16 . EQUILIBRIUM STATE OF A LIQUID 
DROP ON A SOLID SURFACE 
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contact as the intersection of three interfaces — solid-liquid (SL) , 
liquid-vapor (LV) , and solid-vapor (SV) . The balance of these 
surface forces/surface tensions on this three-phase intersection 
is best described by the Young-Dupre equation; (1) 


Y SV " ^SL = Y LV COS 9 (3) 

where 

Y sv = surface tension at 

y gL = surface tension at 

y lv = surface tension at 

0 = contact angle 


When 0 >90°, the liquid is non -wetting and for 90°>e>0°, 
the liquid is wetting but it is nonspreading. When 0 = 0°, the 
liquid wets the solid completely and spontaneously and spreads 
freely over the surface at a rate depending on the liquid viscos- 
ity and solid surface roughness. Nevertheless, it is also possible 
for a liquid to spread and wet a solid with contact angle 0>O° by 
the application of an external pressure or force to the liquid. 

To apply these concepts to adhesive technology, the surface 
and interfacial free energies need to be evaluated. However, there 
are no direct methods available to measure surface free energy, y sv , 
which complicates the analysis. Nonetheless, the contact angle 0 
is a physically measurable parameter and it serves as a convenient 
means for describing the wettability or degree of wetting on solid 
surfaces . 

When a liquid has a zero contact angle with a solid substrate, 
surface tension gradients sometimes exist at the spreading front 
which may either assist or hinder spreading depending upon their 
directions . (2 / 3 ) These surface tension gradients arise from thermal 
gradients, or, in the case of liquids, from the presence of a more 
volatile component of different surface free energy (e.g., a trace 
of impurity) . The effect of the surface tension gradients on the 
spreading rates are relatively small. However, it might have an 
effect on the redistribution of the wetting liquid after it wets the 
solid surface. 

The topography of the substrate surface also has an effect on 
the kinetics of wetting. A liquid can spread along fine pores, 
scratches and other inhomogeneities via capillary action. This 
spontaneous spreading occurs even if the liquid is non-wetting on a 
smooth, planar surface. A fifty per cent increase in spreading rate 
was reported by Bascom et al'2) and Cottington et al.(4) Again, 
topography plays a minor role in affecting the rate of spreading, 
especially if the liquid is forcibly spread over the substrate sur- 
face. Nonetheless, topography does affect the redistribution of 
the wetting liquid after its initial "wet-out." 


solid-vapor interface 
solid-liquid interface 
liquid-vapor interface 
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Much work (5-12) has been done to obtain quantitative ex- 
pressions for wetting upon various types of substrate surface 
topography. Fundamental understanding of the importance of 
wetting kinetics upon wetting processes (e.g., a process to elim- 
inate trapped air at the interface during wetting) and interfacial 
bond strength have yet to be studied. 

Dynamic wetting occurs when the liquid or the solid or both 
are kept in motion relative to the other phase throughout the 
wetting process. This process is best described by a sessile drop 
placed on an inclined plane where two contact angles are present. 
As shown in Figure 17, the advancing angle is the larger, frontier 
angle as the liquid is progressing to "wet" the solid surface. 

The advancing angle reflects the tendency of the liquid to wet a 
virgin surface. The receding angle is the smaller, tail end 
angle where the liquid is retreating from the "wetted" surface. 

The receding angle is a measure of the ability of the solid sur- 
face to remain wetted. 



FIGURE 17. A SESSILE DROP ON AN INCLINED PLANE 

(GRADIENT = tan<f> ) SHOWING THE ADVANCING 

0 = AND RECEDING 0 ANGLES 
a. r 


The dynamic wetting angles best characterize the condition 
where the wetting liquid is forcibly spread onto the substrate 
surface. Moreover, the advancing angle is very sensitive to 
surface contamination and surface heterogeneity. This property 
enables the characterization of the substrate surface and provides 
insight on the wetting process. The receding angle can be related 
to the uniformity of the coating after the liquid has wetted the 
substrate surface. 
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In addition, wetting hysteresis can be obtained by the differ- 
ence between the advancing angles and the receding angles. The 
hysteresis of wetting is due to the presence of a large number of 
metastable states which differ slightly in energy at the surface 
and are separated from each other by small energy barriers . These 
metastable states are generally attributed to either surface rough- 
ness or surface chemical heterogeneity, or both. Surface hetero- 
geneity can also arise from impurities concentrated at the surface, 
from crystal imperfections or from differences in the properties 
of different crystal faces. In general, surfaces are both rough 
and heterogeneous. It is possible to compare surface roughness and 
heterogeneity by measuring the hysteresis of wetting. 

Experimental Techniques 


The surface energetics of fibers can critically affect the 
wetting behavior and the end-use performance of fibrous materials, 
e.g., in composite structures. The measurement of spontaneous sur- 
face wetting is one of the few experimental techniques available 
for the study of the solid-surface energetics. 

Assessing contact angle measurements on filamentous material 
is a very difficult experimental task due to the fiber geometry. 

A preferred technique in direct contact-angle measurement for fibers 
is a derivative of the tilted plate concept. ' 13 ) This involves a 
small reservoir of liquid which is pierced by the fiber and an 
arrangement whereby the fiber can be tilted relative to the liquid 
surface until the meniscus disappears. The contact angle is then 
measured from the fiber to the liquid surfaces. However, considerable 
precaution must be taken to obtain the true contact angle. This ex- 
perimental technique is extremely difficult to apply on fibers with 
small-cross sectional areas. 

An alternate experimental approach which does not require the 
direct measurement of contact angles is the wetting force measure- 
ment based on the Wilhelmy balance principle . (14 ) This technique 
is deemed to be promising with filamentous materials, especially 
with small denier (fine) fibers because of the development of high 
sensitivity microbalances and the ability to assess both equilibrium 
and dynamic wetting data. For fibrous composites, dynamic wetting 
data is very significant because of its ability to prescribe the 
conditions for prepregging processes. 
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The wetting force measurement, according to Wilhelmy, (14) 
is the force exerted on a vertical rod (fiber) inserted into a 
liquid. This force is expressed by: 

F w = * D Y lv cos 0 (4) 

where 

F w = wetting force 

D = diameter of the fiber 

Y = surface tension 
LV 

6 = contact angle 


A mathematical derivation of this expression was performed 
by A. J. G. Allan. (15) The balance of forces acting on the fiber 
before and after contact with liquid is illustrated in Figure 18 . 

An important feature of this technique is its applicability to both 
equilibrium and dynamic conditions. The advancing wetting force 
can be obtained during the fiber penetrating into the wetting liquid. 
The equilibrium wetting force is obtained during the stable condition 
when the fiber is not in motion with respect to the liquid. The re- 
ceding wetting force can be measured during the pulling out of the 
fiber from the wetting liquid. The corresponding wetting force ex- 
pressions are: 


II 

£ 

Pm 

7T 

D 

y lv 

cos 

0 

(5) 

F = 
a 

7T 

D 

y lv 

cos 

0 

a 

(6) 

F = 
r 

7T 

D 

y lv 

cos 

0 

r 

(7) 


where 

F w and 0 are the equilibrium wetting force 
and contact angle respectively 

F and 0 are the advancing wetting force 
a a and contact angle respectively 

F and 0 are the receding wetting force 
r r and contact angle respectively 

By measuring the wetting force, F^, and the surface tension 
of the wetting liquid, y , work of adhesion, W a can be determined. 

LV 


or 


W A = 


LV 


+ Y lv COS 


( 8 ) 


W, 


LV 
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(9) 
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(a) Filament suspended 
freely 


F b M f g 

(b) Filament partially 
immersed 


where : 


F b = buoyancy force Y LV = surface tension of 

M^g = gravitational force wetting liquid 

F^ = wetting force 0 = contact angle 

D = diameter of fiber 

assuming F^ is negligible 

then 


F 


W 



TT D Y COS 0 
ij V 


FIGURE 18. FORCES ACTING ON A SUSPENDED FILAMENT BEFORE AND 
AFTER CONTACT WITH WETTING LIQUID 


44 



Furthermore, by attaching a microprocessor to the micro- 
balance and controlling the fiber immersion/emersion operation, 
a continuous monitoring of the advancing wetting force and receding 
wetting force along the fiber can be achieved. This would enable 
the determination of fiber surface heterogeneity and provide in- 
sight into the uniformity of the liquid film adhered to the fiber 
surface. Moreover, this information also shows the frequency of 
occurence of critical surface contamination which prevents wetting/ 
bonding. 

Critical surface contamination is present in those areas of 
the surface that has a surface energy lower than the wetting liquid. 
This lower surface energetics create a thermodynamically unstable 
condition for the higher surface energy liquid to "wet" the solid 
surface. Therefore, a localized rejection of the wetting liquid 
occurs and this leads to the formation of voids during composite 
fabrication . 

In addition, the hysteresis of wetting can also be studied by 
this technique. A plot of the immersion f orce/emersion force along 
the fiber immersion length can be used to map the wetting hysteresis. 
Figure 19 shows a typical plot and identifies the hysteresis of 
wetting . 



FIGURE 19. HYSTERESIS OF WETTING MEASURED BY PLOTTING ADVANCING 
AND RECEDING WETTING FORCE VERSUS IMMERSION 
FIBER LENGTH 
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Although the Wilhelmy wetting force measurement can provide 
much information, one experimental parameter can affect the accuracy 
of all these data. (16) The rate of the fiber penetrating into and 
receding out of the wetting liquid critically affects the measur- 
ables. As shown in Figure 20, a 90° contact angle (total disappear- 
ance of the meniscus) is created by this dynamic phenomena at the 
"critical penetration velocity." This phenomenon cannot be totally 
eliminated. However, the degree of meniscus depression/elevation can 
be minimized and controlled by operating at a much slower fiber 
immersion/emersion rate than the critical velocity and maintaining 
this rate throughout all studies. This technique has been used by 
many authors (17-23 ) an £ information on small denier fibers, such as, 
carbon fibers and polyester fibers, have been reported. 


v<v v=v v>v 

C C c 

I i I 



where 


V = fiber penetration velocity 

V = critical velocity at which meniscus 
disappears as in (c) 


FIGURE 20. EFFECT OF FIBER PENETRATION VELOCITY ON CONTACT ANGLES 


Carbon Fiber Surface Energetics 


Hammer and Drzal(19) studied the surface energetics of Hercules 
Type A carbon fibers via the micro-Wilhelmy wetting force technique. 
Correlations were found between the oxygen concentration measured by 
x-ray photoelectron spectroscopy (XPS) and the change in polar/dis- 
persive ratios of the carbon fibers. This indicated that surface 
treatments of carbon fibers which promote better fiber matrix adhesion 
resulted from an increase in the surface free energy. This increase 
in surface free energy was due to the addition of surface oxygen con- 
centration during the surface treatment process . Some of the contact 
angle data are shown in Table III while Table IV lists the measured 
surface free energies of carbon fibers with their polar and dispersive 
components . 



TABLE III. CONTACT ANGLES FOR CARBON FIBERS* 


Liquid 


T LV 

(mJ/m 2 ) 


AU 


AS 


As Received 300°C V.T. As Received 300°C V.T. 


Water 

72.8 

0 

H 

± 

4° 

44° 

± 

4° 

29° 

± 

4° 

33° 

± 

4 ° 

Glycerol 

64.0 

33° 

± 

2° 

40° 

± 

1° 

25° 

± 

1° 

42° 

± 

9° 

Formamide 

58.3 

37° 



24° 

+ 

7° 

22° 

± 

5° 

27 ° 

+ 

3° 

Methylene Iodide 

50.8 

23° 

± 

5° 

41° 

± 

11° 

23° 

± 

5° 

30° 

± 

12 

Ethylene Glycol 

48.3 

35° 

± 

6° 

30° 

± 

6° 

31° 

± 

7° 

29 ° 

± 

2° 

Br omonaphtha lene 

44.6 

29° 

± 

3° 

20° 

± 

4° 

26° 

+ 

3° 

21 ° 

± 

3° 

Polypropylene 

31.3 

3° 

± 

3° 

25° 

± 

1° 

3° 

± 

3° 

25 ° 

± 

1° 

Glycol PG 1200 














n-Hexadecane 

27.6 

18° 

± 

8° 

27° 

± 

7° 

22° 

± 

7° 

23 ° 

± 

7° 


Y lv = surface tension of liquid 

AU = Hercules Type A high strength carbon 
fiber non-surface treated 

AS = Hercules Type A high strength carbon 
fiber surface treated 

V.T. = vacuum treated 


*From Hammer and Drzal (ref. 19) . 
Reproduced by permission. 
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TABLE IV 


SURFACE FREE ENERGIES OF CARBON FIBER* 


Fibers 


Surface Free Energies 

As Received 

P 

Y sv (mJ/m 2 ) 

Y g v (mJ/m 2 ) 

Y gv (mJ/m 2 ) 

AU 

23.6 ± 2.6 

27.4 ± 0.3 

51.0 ± 2.6 

AS 

30.0 ± 1.7 

26.4 + 0.1 

56.4 + 1.7 

300 °C V.T. 

AU 

24.1 ± 1.4 

26.3 ± 1.1 

50.4 ± 1.4 

AS 

26.8 + 1.4 

26.0 ± 1.2 

52.8 ± 1.4 


p 

Y gv - Keesom-polar surface free energy of solid. 

Y^ v - London -dispersion surface free energy of solid. 

p D 

Y gv - Surface free energy of solid, y sv = Y sv + Y sv 
AU - Hercules Type A carbon fiber — non-surface treated. 
AS - Hercules Type A carbon fiber — surface treated. 

V.T. - Vacuum treated. 


*From Hammer and Drzal (ref. 19). 
Reproduced by permission. 
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Kaelble et a l ( 21 ' 2 ^" 2 5 ) have intensively studied carbon fiber 
surface energetics via the micro-Wilhelmy balance technique. Some 
of his data on the contact angles of carbon fibers are shown in 
Table V. Kaelble also extended the information on the Keesom- 
polar and London-dispersive surface free energy with the modified 
Griffith fracture energy. This enabled the prediction of the inter- 
facial conditions to promote good bond strength and environmental 
stability. This expression was done through the extension of the 
adsorption theory of interfacial interactions . ( 2< ^ ' The critical 
value of the Griffith fracture stress, <? c under normal plane stress 
loading can be expressed by the following relationships: 


a 


c 




2 2 j, 

(R -R P > 0 
o 


( 10 ) 



H = 
K = 



a 2 and S 2 ~ 
For air: 

For water: 


N> 

1 

0 m 


(ii) 

0.25C(a 1 - a 3 ) 2 

+ (b 1 - e 3 ) 2 : 

(12) 

(« 2 - H) 2 + (B 2 

- K) 2 

(13) 

0.5 (a ^ + a 3 ) 


(14) 

0.5 (3 1 + B 3 ) 


(15) 

(Y d ) ^ 
VT LV ; 


(16) 

(Y P ) % 
V T LV' 


(17) 

(Y d ) H 
T sv' 


(18) 

(y p 

w sv 7 


(19) 


surface properties of the immersion phase 
at the crack tip 

a 2 = 6 2 = 0 

a 2 = 4.67 (mJ/m 2 ) 32 

B 2 = 7.14 (mJ/m 2 ) 35 
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TABLE V. WILHELMY BALANCE MEASUREMENTS OF HERCULES 
TYPE AU AND TYPE AS CARBON FIBER AT 23 °C* 


AU AS 


Liquid 

y lv 

(mJ/in 2 ) 

F 

(yg) 

COS0 

a 

0 

a 

W a 

(mJ/m 2 ) 

Fa 

(yg) 

COS0 

a 

0 

a 

w 

a 

(mJ/in 2 ) 

Water 

72.8 

181 

0.8299 

33.9° 

133.22 

172 

0.8966 

26.3° 

138.07 

Glycerol 

64.0 

182 

0.9492 

18.3° 

124.75 

161 

0.9547 

17.3° 

125.10 

Formamide 

58.3 

131 

0.7500 

41.4° 

102.03 

120 

0.7811 

38.6° 

103.84 

Di i odome t h an e 

50.8 

115 

0.7556 

40.9° 

89.18 

102 

0.7620 

40.4° 

89.51 

Ethylene Glycol 

48.3 

137 

0.9467 

18.8° 

94.03 

120 

0.9429 

19.5° 

93.84 

1-Bromonapthalene 

44.6 

126 

0.9430 

19.4° 

86.66 

111 

0.9445 

19.2° 

86.72 

Glycol PG-E-1200 

43.5 

lBO^l.OOO 

0° 

— 

114 

1.000 

0° 

— 

Glycol PG-15-200 

36.6 

106 

0.9667 

14.8° 

71.98 

94 

0.9747 

12.9° 

72.27 

Glycol PG-P-1200 

31.3 

94 (1) 1.000 

0° 

— 

83 (1) 1.000 

0° 

— 

n-Dodecane 

25.4 

71 

0.933 

21.1° 

49.10 

63 

0.9413 

19.7° 

49.31 

Ethanol (abs) 

22.4 

— 

— 

— 

— 

56.3 

0.9413 

19.7° 

43.49 


Y - surface tension of liquid 
LV 

F - advancing wetting force 
a 

0^ - advancing contact angle 

W - work of adhesion 
a 

AU - Hercules type A non-surface treated carbon fiber 
AS - Hercules type A surface treated carbon fiber 

(1) Wetting conditions with cosine 0 = 1.0. 

*From Kaelble et al (ref 21) 

Reproduced by permission. 
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Through this thermodynamic analysis, Kaelble predicted that 
in order to have a moisture insensitive interface, a non -polar sur- 
face for the fiber and the matrix C3i, 3 3 ^1 . 0 (mJ/m 2 ) hi was neces- 
sary. To assure spontaneous wetting to enhance interfacial bond- 
ing, an increase in London -dispersive contribution to the fiber 
surface free energy C« 3 ^6 . 0 (mJ/m 2 ) is required* 

Penn et al( 22 ) studied the advancing and receding contact 
angles of carbon fibers and the effect of wetting hysteresis. 

The advancing angle is used to demonstrate how readily the liquid 
wets the fiber surface; the receding angle is to illustrate 
how well the fiber surface remains wetted. The data is presented 
in Table VI. According to Penn, hysteresis of wetting would 
disappear below a critical surface tension of the wetting liquid. 

In addition, Penn also tried to correlate the wetting properties 
of carbon fibers to composite shear strength . ( 22 ) Penn sized the 
carbon fiber with a silicone fluid to create a hydrophobic surface 
and also had a hydrophilic size to compare the effect of sizings 
to composite properties both unconditioned and conditioned. 

Some of these data are shown in Table VII. However, the re- 
ported composite data do not correlate with the wetting angles 
measured. The hydrophilic sized fiber showed the smallest con- 
tact angles (best wetting) but it had the lowest short beam shear 
strength both conditioned and unconditioned. Furthermore, silicone 
fluid has a poor compatability with epoxy and during prepregging, 
the silicone fluid could diffuse away from the fiber surface to 
the prepreg surface. This effect would be more pronounced with 
filament winding procedures. This is because the filament winding 
resin systems are more fluid and molecule mobility is enhanced. 

The similarity of the composite short beam shear strength of the 
silicone fluid sized fiber and the unsized fiber could be the 
reflection of the disappearance of the silicone size from the fiber 
surface. Moreover, Penn stated that silicone fluid was found on 
the surface of the laminates.* 

Overall, much work has been done in measuring contact angles 
of carbon fibers and theoretical analyses have been made. However, 
reported values are generally scattered and the conditions at 
which these measurements were made are not generally clear. Above 
all, fiber origins, types and properties have not been generally 
characterized and specified. Therefore, attempting to correlate 
the reported data is difficult and unclear. This report is in- 
tended to clarify some of the gray areas and to try to correlate 
wetting properties of carbon fibers to fiber matrix adhesion and 
composite properties. 


*Private communication with Dr. Lynn Penn . 


52 



TABLE VI. ADVANCING AND RECEDING CONTACT 
ANGLES OF CARBON FIBERS* 


Y Carbon Fiber ^ 

LiV — — 


Liquid 

(mJ/m 2 ) 

COS0 _ 

a 

!a 

COS0 r 

fr 

Water 

72.4 

-0.067 

93.8° 

0.713 

44.5° 

Formamide 

55.6 

0.445 

62.9° 

0.562 

55.8° 

Methylene Iodide 

48.6 

0.824 

34.5° 

0.917 

23.5° 

Ethylene Glycol 

47.5 

0.538 

57.5° 

0.895 

26.5° 

Br -Naphthalene 

45.3 

0.957 

16.9° 

0.957 

16.9° 

Dimethyl Acetamide 

36.7 

0.993 

6.78° 

0.993 

6.78 

Hexadecane 

28.4 

0.964 

15.4° 

0.964 

15.4° 


(1) Origin and type of carbon fiber unknown. 
Y - surface tension of liquid. 

LV 

6 - advancing contact angle. 

0 r - receding contact angle. 


*From Penn and Miller (Ref. 22) 
Reproduced by permission. 
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TABLE VII 


EFFECT OF WETTABILITY ON CARBON FIBER/EPOXY COMPOSITE SHEAR STRENGTH* 


Short Beam Shear (AS TO D2344) 


14 Days Water Boil 


Fiber Variants 

0 

a 

0 

r 


RT 82 . 2 °C (180 °F) 



RT 82.2° 

C (180°F) 





MPa 

ksi 

MPa 

ksi 

MPa 

ksi 

MPa 

k"si 

AS- 4 

(unsized) 

52.3° 

22.0° 

98.614.8 

14.310.7 

78.612.1 

11.4+0.3 

80+2.1 

11. 610. 3 

58.610.7 

8.510.1 

AS -4 

(hydrophobic size) 
(silicone fluid) 

78.6° 

24.6° 

90.3+2.8 

13.1+0.4 

76.6+1.4 

11.1+0.2 

81.412.8 

11.810.4 

58.6+2.1 

8.510.3 

AS -4 

26.7° 

0.00° 

69.710.3 

10.110.05 57.2+4.1 

8.310.6 

42.8+4.1 

6.210.6 

33.8+2.1 

4.910.3 


*From Desalvo et al (ref 23) 
Reporduced by permission. 
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STRESS TRANSFER FROM MATRIX TO FIBER 


When a tensile load is applied to an elastic matrix, the 
matrix will elongate to a value proportional to the load. How- 
ever, if fiber with a higher elastic modulus is embedded in the 
matrix, the fiber will restrict the matrix from freely elongating 
in the region of the fiber. Shear stresses are produced on planes 
parallel to the axis of the fibers in the direction of this axis. 
The tensile stresses in the fiber are introduced by these inter- 
facial shear stresses. These stress patterns can be identified 
and evaluated by birefringence techniques.^' 2 ) 


A number of models have been proposed for calculating the 
longitudinal tensile stress distribution in an embedded fiber. 

The models of Cox, (3) Dow( 4 ) and Rosen( 5 ) considered the case of 
an elastic fiber in an elastic matrix. If the matrix deforms elas- 
tically, the interfacial shear stresses peak at a maximum value near 
the fiber ends and then decay rapidly toward the center length of 
the fiber. Conversely, the tensile stresses in the fiber will build 
to a maximum value provided that the fiber is sufficiently long. 

The tensile strain of the fiber and the tensile strain of the matrix 
are nearly equal in the mid-length. 


When the load is applied to the fiber by shear stress at the 
interface, physically, at equilibrium, one must always have; 


dP 

dx 


IT 


d f T x,r= J 3d f 


( 20 ) 


where 

x = distance from end of fiber 
P = tensile load in the fiber 
d = diameter of the fiber 

= the shear stress evaluated 
x,r=*sd f at the fiber surface 

Generally, r in the matrix is a function of both x and r. 

The magnitude of the interfacial shear stress is largely dependent 
on the fiber content, on the magnitude of the applied load, and on 
the physical properties of the filament, matrix and interfacial 
bonds . 

In order to obtain an expression for the shear stress, the 
following assumptions are made. (3-5) 
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1. Fiber and matrix both behave elastically. 

2. Fiber is surrounded by a cylinder of matrix material. 

3. The bond between fiber and matrix is perfect throughout 
any deformation. 

4. Straight radial lines in the fiber and matrix remain 
straight lines after deformation. 

The shear stress between fiber and matrix could then be ex- 
pressed as follows: (3-5) 


T 


x 



( 21 ) 


where : 



( 22 ) 


and 

r = interfacial shear stress between fiber and matrix 
x in psi 

A = cross-sectional area in in 2 
E = Young's modulus in psi 
G = shear modulus in psi 

X = distance from end of fiber, in inches 
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d^ = diameter of fiber, in inches 

L = effective length for load differential between 
fiber and matrix 



o a = applied field stress in psi 

Subscripts : 

f = fiber 
m = matrix 

Hence, the tensile load in the fiber at distance X is given by: (3-5) 

p f = P f ,x=o + * d f / T x ** < 24 > 

Therefore, the tensile stress in the fiber at position X is 
obtained by integrating as follows : (3-5) 
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where : 



The transfer length required to build up a maximum stress 
equal to a fraction 0 of the stress in an infinitely long fiber 
can also be calculated by the equation: (3-5) 


^ = A cosh - 1 ( -±- ) 
d f A \ 1-flJ / 


(27) 


and the load transfer pattern is demonstrated in Figure 22. 


Matrix 


Matrix 



L, 


Reference 
Line o 1 




A A 


■H 


f Hh Hh f '° 

L L r. 

C C 

a. Elastic Fiber 
Elastic Matrix 



b. Elastic Fiber 
Plastic Matrix 


FIGURE 22. INTERFACIAL SHEAR STRESS PATTERNS FOR THE 
ELASTIC AND PLASTIC MATRIX CASES 
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Prom equation (21), one can see that a maximum shearing stress 
occurs at the ends of the fiber. When a high tensile load is 
applied, the shear stress is likely to exceed the yield stress of 
the matrix. If the matrix is unable to plastically deform, either 
the matrix fails or an interfacial slip occurs between the fiber 
and the matrix. If the matrix is able to deform, then the inter - 
facial shear stress never rises above the yield shear strength of 
the matrix and equation (21) becomes irrelevant. Therefore, the 
plastic flow model of Kelly and Tyson (6) appears to be more approp- 
riate . 


In the case of a plastic matrix, a yield shear strength x y is 
constant throughout the deformation. Then one can obtain a tensile 
stress distribution in the fiber by integrating equation (20) : (6) 


P = a 


■ / 


ird-t dx 

f y 


(28) 


= Trd,-T x 

f y 


Equation (28) shows that the tensile load builds up linearly 
and a transfer length of the fiber is obtained. (.6) 


L 


c 


d 


f 


4t 

y 


(29) 


where 

CT , 


or 

T 

y 

= tensile strength of the fiber at transfer length 
= fiber diameter 
= inter facial shear strength 

= critical length of the fiber required for load transfer 
from a fiber end 


a d 
f f 


2 Z 


(30) 


Z 

c 


critical fiber length required for load transfer 
Z = 2 l 


c c 


/ 


and the load transfer mechanism for the plastic deformable matrix 
is shown in Figure 22 (b) . 


61 


References to Stress Transfer 


(1) L. B. Greszczuk, "Micromechanics Failure Criteria for Composites 

Subjected to Transverse Noarmal Loadking," AIAA Journal, 9, 

No. 7 (1971), 1274. 

(2) A. Voloshin, "Stress Field Evaluation in Photoelastic Anisotropic 

Materials: Experimental Numerical Technique," J. Comp. Mat'l., 
14 (1980), 342. 

(3) H. L. Cox, "The Elasticity and Strength of Paper and Other 

Fibrous Materials," Brit. J. Appl. Phys . , 3 (1952), 72. 

(4) N. F. Dow, "Study of Stresses Near a Discontinuity in a Filament- 

Reinforced Composite Metal," General Electric Report TIS 
R635D61, AD-414673 (1963). 

(5) B. W. Rosen, "Tensile Failure of Fibrous Composites," AIAA J., 

2, No. 11 (1964), 1985. 

(6) A. Kelly and W. R. Tyson, "Tensile Properties of Fibre-Reinforced 

Metals: Copper /Tungsten and Copper/Molybdenum, " J. Mech. 

Phys. Solids, 13 (1965), 329. 


62 



METHODS FOR MEASURING INTERFACIAL BOND STRENGTH 


Fiber Pullout Test 

Broutman^^ developed the single filament fiber pullout test 
to measure the interfacial joint strength and interfacial frictional 
strength. A typical load -displacement curve resulting from such 
a test is shown in Figure 23. The geometry of the specimen for the 
fiber pullout test can either be a but ton -type sample as shown in 
Figure 24, or a fiber end of controlled length embedded into a 
bulk matrix as in Figure 25. The relationship which assumes the 
shear stress to be uniformly distributed along the interface is: 


T 


'm 


2 irr £ 


or 

m 


2 £ 


(31) 


where 

t = average shear strength of joint 
= maximum load applied to fiber 
r = radius of fiber 
£ = embedded fiber length 

cr = maximum stress applied to fiber 


The embedded fiber length is influenced by the fiber strength 
so that the maximum embedded length to be used is determined by: 


where 


£ 

max 


a 


ult 
2 T 


r 


(32) 


a lt = ultimated tensile strength of the fiber 

. ( 2 ) 
The theoretical analysis of this test was performed by Lawrence ' 

and Barto^(3) w ho concluded that this test was sensitive to interfacial 

adhesion properties. 

However, this test is not applicable to carbon fiber/epoxy 
joint strength analysis because of the small diameter of the carbon 
fiber and the relatively high interfacial bond strength. Using the 


63 




FIGURE 23. A TYPICAL LOAD -DISPLACEMENT CURVE 
FOR FIBER PULLOUT TEST* 


*From Broutman (ref. 1) 
Reproduced by permission. 
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FIGURE 


FIGURE 
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Matrix 
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24. BUTTON-TYPE SPECIMEN FOR FIBER PULLOUT TEST 



25. BLOCK-TYPE SPECIMEN FOR FIBER PULLOUT TEST 
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based on P and % as in equation (33) . 

HI/ C C 


P 

_ m/C 

7rd& 

C 


( 33 ) 


where : 

P m,c = fiber pullout force at maximum embedded length 

= maximum embedded length which allows fiber pullout 


This relation was not determined by Favre possibly because of the 
difficulties in controlling and fabricating ultrathin resin lozen- 
ges. Some of his data is shown in Tables IX and X. 


TABLE IX. VALUES OF THE ADHESION STRENGTH a 
FOR AS -RECEIVED FIBERS* A 


Fibre 

Type 

Modulus (GNm -1 ) 

Matrix • 

o a Nmm~* 

Counauld5 

HT 

270 

Polyester A 

5.4 


(long staple) 


Polyester B 

10.4 




Epoxy 

9.7 




Polyimide 

17.1 


HTS 

283 

Epoxy 

57.3 


(continuous tow) 

HM 

330 (1st batch) 

Epoxy 

5.5 


(continuous tow) 


Polyimide 

12.8 



357 (2nd batch) 

Epoxy 

1.8 


HMS 

(continuous tow) 

350 

Epoxy 

31.4 

Morgan 

11 treated (sample) 

240 

Epoxy 

41.9 

Carbone-Lorraine 

AG (short staple) 

440 

Epoxy 

1.0 


AGT (continuous tow) 

258 

Epoxy 

41.6 

Columbia 

CMR (sample) 

42t 

Epoxy 

8.3 


‘Polyester A Rhodester 1108 (RhSne-Poulenc) 
Polyester B Stratyl A. 30 (St Gobain) 

Epoxy Araldite LY 556 + HT 972 (CIBA) 

Polyimide Kerimid 601 ( Rh3ne -Poulenc ) 

fManuf acturer ' s data 

*From Favre and Perrin (Ref. 4) 

Reproduced by permission. 
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In Table IX, 


F 

0 A ird£ 


where 


F = the force required for fiber pullout 
d = diameter of the fiber 
a = thickness of the pastille 


TABLE X. ADHESION STRENGTH OF TREATED FIBRES* 
(Matrix: epoxy LY 556 + HT 972) 


Fibre Treatment Modulus Tensile strength Nmnr* 

(GNrrr*) (GNm J ) 


Courtaulds HT 

— 

270 

2.56 

9.7 

(long staple) 

benzene washing 

270 

2.56 

7.8 


oxidation (nitric acid) 

— 

2.14 

*2.5 


oxidation (hypochlorite) 

273.5 

2.46 

40.0 


oxidation (hot air) 

— 

1.96 

40.5 


formation of silicon 
carbide 



18.0 

Courtaulds HM 



357 

1.73 

1.8 

(continuous) 

oxidation (hypochlorite) 

270 to 325 

1.00 to 1.80 

up to 26 


formation of silicon 
carbide 

281.5 

0.88 

(70) 


reduction by wet hydrogen 259 
pyrolytic carbon coating 

1.04 

3.8 


•toluene 

340 

1.83 

8.0 

From Favre and 
Reproduced by 

-propylene 

Perrin (Ref. 4) 
permission. 


1.67 

8.0 to 46 
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Single Filament Shear Fracture Energy Test 

Outwater and Murphy (5-7) developed a new test method to meas- 
ure the debonding fracture energy between epoxy resin and single 
filament glass fibers. The authors considered an infinitely long 
fiber of diameter "df" and cross-sectional area "Af" embedded 
in a semi -inf inite solid as in Figure 27. Let P be the tensile 
load at the end of the fiber protruding a distance L from the sur- 
face of the' resin and let the fiber be debonded a distance "x" 
into the matrix and be the energy required to debond the fiber. 

The strain energy, £2 f in the filament is: 

P 2 L /*x (P— TTfd-x) 2 

fl £ = + / dx 

2AE f J 0 2AE f 


where 


x = interfacial shear strength 
= Young's modulus of the fiber 



FIGURE 27. EMBEDDED SINGLE FIBER IN A SEMI-INFINITE 
SOLID UNDER TENSILE LOAD* 

*From Outwater and Murphy (Ref. 5) . 

Reproduced by permission. 
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By equating the energy required to debond an incremental 
length, dx, Gjjirdfdx and the strain energy released by an incre- 
ment dx of debonding. 


(P-xird-x) 2 

- dx 

2AE f 


where A = cross-sectional area, the following equation is 
obtained: 


G 


TT d 

II f 


(P-T1Td f X) 2 

2AE f 


(35) 


Substituting for P = a*A, then 


where 



(36) 


a = stress at the fiber end not embedded in resin 


Directly pulling the fiber makes this measurement experimen- 
tally difficult. A modified technique is to embed the fiber in a 
resin block, severing the fiber at the center of the block and 
then loading the specimen in compression and observing the stress 
required to debond the fiber. A diagram of this method is shown 
in Figure 28. The G^ value was determined from: 



(37) 


where 


= applied stress 
x = debonded length 

and, for initial debonding, where x =0, then 


G 


II 


e 2 e 
r 


f a f 


8 


( 38 ) 


where 


e 


r 
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= compressive strain in the resin 




FIGURE 28 . SPECIMEN FOR MEASURING THE DEBONDING ENERGY 
OF THE FIBER FROM THE SURROUNDING MATRIX* 

*From Outwater and Murphy (Ref. 7) . 

Reproduced by permission. 
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The data obtained from this test for glass fibers in an epoxy- 
matrix was conclusive. However, due to the small diameter and the 
brittle nature of the carbon fiber, this test is not applicable 
for a study of carbon fiber /epoxy adhesion properties. This test 
requires the observation of initial debonding during compression 
of the specimen and the small diameter of a carbon fiber makes the 
observation difficult if not impossible. 

Single Filament Interfacial Shear and Transverse Tensile Debonding 
Strength Determination " 


Broutman (1 / 8 , 9) has also developed tests to characterize the 
interfacial shear debonding strength and interfacial transverse 
tensile strength. Interfacial shear strength is determined by 
embedding a single filament in the center of an epoxy casting as 
shown in Figure 29. Upon subjection to compression, the fiber 


.27 cm — *| 



T 


3 . 81cm 




FIGURE 29. INTERFACIAL SHEAR STRENGTH SPECIMEN 


will debond at the interface by shear. The shear stress can be 
determined from: 


x = 2 . 5 ct 


(39) 


where 


t = interfacial shear stress 
o = axial compression stress in specimen 
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Interfacial transverse tensile strength was determined in a 
similar manner. The specimen for this test is shown in Figure 30. 



FIGURE 30. INTERFACIAL TRANSVERSE TENSILE STRENGTH SPECIMEN 


Due to the curvature of the matrix and the difference in the 
Poisson's ratios of the fiber and the matrix, the fiber will 
debond at the interface through transverse tensile stress. This 
transverse tensile stress can be determined from: 

S = - < V (|J m~ ti f )E f (40) 

< 1 + V E f + 

where 

S = transverse tensile stress 
a m = axial compression stress on specimen 
y = Poisson's ratio 
E = elastic modulus 

Subscripts: f = fiber 

m = matrix 

These tests have been used to evaluate carbon fiber/epoxy ad- 
hesion properties by Hawthorne ) an( j Mullin. (H) However, they 
both concluded that these tests were not applicable to carbon fibers 
due to the low compression strain of the carbon fibers. Compression 
fracture and local buckling occured prior to the debonding of the 
fiber . 
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Single Filament Critical Length De termination 

The single filament critical length determination is an experi- 
mental technique derived from the load transfer analysis discussed 
under "Stress Transfer from Matrix to Fiber." Kelly (16) illustrated 
the multiple fracture phenomenon in a system consisting of tungsten 
wires embedded in a copper matrix. After deformation, the copper 
matrix was dissolved and the wire fragments extracted. On analysis 
of the fragment length distribution, it was found that its bounds 
were approximately a c and £ c /2; where j l c was the critical fiber 
length and was determined froio: 


i = — 
c 2 T 

where 

Of = tensile strength of fiber at Z c gauge length 
d f = diameter of fiber 
r y = interfacial shear strength 

Equation (41) is a rearrangement of equation (30) . Recalling the 
definition of critical fiber length, this result is not unexpected. 
The longest fibers capable of surviving in this load situation will 
have a length just under a ; since longer fibers, by definition, 
must fracture. c 

Analogous experiments were run by Ongchin et al, (18) McGarry(19) 
and Fraser et al(12) for glass fibers in both thermoplastic and 
thermoset (epoxy) matrixes. Conclusive results were obtained. 
Wardsworth and Spilling, (17 ) ishikawa et al(13) and Drzal et al (14*15) 
applied this technique to carbon fiber in an epoxy matrix. 

Drzal (14) embedded a single filament of carbon fiber into an 
ASTM two-and-one-half inch length epoxy casting as in Figure 31 . 

These specimens were then subjected to a tensile deformation. The 
schematic of the fiber fracture and stress distribution is presented 
in Figure 32 (taken from Ref. 14). Figures 33 and 34 show the 
transmitted polarized light micrographs of a non-surface treated 
(fiber A) and a surface treated (fiber B) fiber during tensile de- 
formation and their corresponding shear strength data are presented 
in Table XI. The effect of carbon fiber surface treatments on ad- 
hesion is shown clearly by the histogram of the aspect ratios of 
the two types of fibers in Figure 35. 

In addition, the failure mechanisms during debonding can also 
be observed and characterized. (15) The epoxy sized carbon fiber 
promoted a higher interfacial bond strength. The higher joint 
strength created matrix cracking after fiber fracture whereas no 
matrix cracking was observed with the non-surface treated, unsized 
fiber and the surface treated, unsized fiber. 
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□ 

Section A-A 


FIGURE 31. SCHEMATIC DIAGRAM OF SINGLE FIBER INTERFACIAL 
SHEAR STRENGTH SPECIMEN* 


*From Drzal et al (Ref. 14) . 
Reproduced by permission. 





a. At low initial stress values 

b. After fiber flaws and defects have broken. 

c. After critical length has been reached. 

FIGURE 32. FIBER FRACTURE AND STRESS DISTRIBUTION* 

*From Drzal et al (Ref. 14) . 

Reproduced by permission 
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FIBER A / EPON 828 - mPDA 



FIGURE 33. TRANSMITTED POLARIZED LIGHT MICROGRAPH OF A 
TYPICAL FIBER A FRACTURE AS A FUNCTION OF 
INCREASING STRAIN* 


*From Drzal et al (Ref. 14). 
Reproduced by permission. 




FIGURE 34. TRANSMITTED POLARIZED LIGHT MICROGRAPH OF A 
TYPICAL FIBER B FRACTURE AS A FUNCTION OF 
INCREASING STRAIN*. 


*From Drzal et al (Ref. 14) . 
Reproduced by permission. 
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TABLE XI. SINGLE FILAMENT INTERFACIAL SHEAR 

STRENGTH 


OF FIBER A U ' 

AND FIBER B 12 '* 




Fiber A 

(1) 

Fiber B 

(2) 

Total Specimens 

15 


23 


Total Fragments 

669 


485 


Weibull a (eqn. 3) 

3.074 


3 .308 


Weibull 3 (eqn. 4) 

48.50 


109.05 


Fiber Tensile Strength 

2959 

MPa 

2614 

MPa 

(1" gage length) 

(429,000 

psi) 

(379,000 

psi) 

E828/mPDA Yield Strength 

40 

MPa 

40 

MPa 


(5,800 

psi) 

(5,800 

psi) 

Interfacial Shear Strength (t) 

40 . 86±26 . 21 

MPa 

15 . 6±8 . 89 

MPa 

(eqn. 1) 

(5, 925±3800 

psi) 

(2 , 262±1290 

psi) 


(1) As— 1 — surface treated carbon fiber, unsized. 

(2) AU-1 — virgin carbon fiber with no surface treatment and unsized. 

*From Drzal et al (Ref 14) . 

Reproduced by permission. 
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Total number of breaks as a function of length 
to diameter ratio. 


FIGURE 35. HISTOGRAMS OF ASPECT RATIOS FOR FIBER A AND B* 


*From Drzal et al (Ref. 14). 
Reproduced by permission. 






Overall, this test has shown some conclusive results and pro- 
vided much information on the characterization of fiber/matrix ad- 
hesion. This technique has been adapted by the author to further 
study the effect of interphase properties on carbon fibers epoxy ad- 
hesion. 


Microdebonding Test 

Mandell and McGarry(20) developed a microdebonding test for 
in-situ measurement of fiber/matrix bond strength in fiber com- 
posites. The test involved the compressive loading of a fiber on 
a polished specimen surface to produce debonding as shown in Fig- 
ure 36. Glass, aramid and carbon fibers in epoxy composites have 
been tested. Initial results of these evaluations, especially 
glass, on unconditioned and moisture conditioned samples were 
promising. For carbon fiber, experimental difficulties were encoun- 
tered due to the skin-core properties of the fiber. Refinements to 
simplify interpretation of the debonding force in terms of inter- 
facial shear strength and improvement on the apparatus were required 
in order to make the test more reproducible for carbon fibers. 



FIGURE 36. MICRODEBONDING TEST* 

*From Mandell and McGarry (Ref. 20) 

Reproduced by permission. 
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EXPERIMENTAL 


INTRODUCTION 

After reviewing the information available in the literature, 
it is apparent that little effort has been spent on the correlation 
of theoretical analyses with experimentally determined composite 
properties; and there is limited evidence that the nature of the 
interphase can be used to optimize composite properties. In order 
to clarify some of these uncertain areas, the following experiments 
were designed to examine the relationship between wetting and in^ 
terfacial shear strength. Also, various finish variants were 
formulated to demonstrate the effect of interphase properties on 
composite performance. 


MATRIX RESIN 

For carbon fiber composites, diglycidyl ether of bisphenol A 
(DGEBA) epoxies are commonly used in the matrix resin. In order 
to have a typical matrix that is suitable for both critical 
length determination and composite fabrication, system "F" was 
formulated for the evaluations reported herein. This resin is 
prepared by mixing 100 parts of Epi-Res 508* and 22.5 parts of 
Epi-Pure 841.* A unique feature of this system is its high ten- 
sile elongation-to-break value of 8%. A cure schedule of two 
hours at 93°C was used combined with a post-cure at 150 °C for two 
hours. The tensile properties of system "F" castings, evaluated 
at different cure schedules, are shown in Table XII. This high 
elongation-to-break property of system "F" allows multiple frac- 
ture of the embedded fiber without matrix failure which enables 
critical length determination. Some of the other properties ob- 
tained from castings of system "F" are presented in Table XIII. 

Finish Variants 

To study the effect of fiber/matrix interphase on composite 
performance, four finish variants were developed and applied to 
Celion® 6000 (C-6K) unsized fiber. These four variants are: 

1. C-6K/SR - RTV 615**silicone rubber 

2. C-6K/FRE 25 -An organophosphazene***elastomer blended with 

Epon 828 at 25 parts elastomer to 100 parts 
epoxy 


* Product of 
** Product of 
***Product of 


Celanese Plastics and Specialties Co., 
General Electric Co. 

Firestone Co. 


Louisville, Ky. 
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TABLE XII. TENSILE PROPERTIES OF SYSTEM "F" AT 
VARIOUS CURE SCHEDULES 


Composition: 


EPI-RES 508* 100 

EPI-CURE 841* 22.5 


Cure Post 


Schedule Cure 


Sample 

Time 

ThtT 

Temp 

■Pci 

Time 

ThtT 

o CD 

0 3 

1 

2 

85 

2 

150 

2 

2 

85 

3 

150 

3 

2 

85 

4 

150 

4 

2 

85 

2 

160 

5 

2 

85 

2 

170 

6 

2 

85 

2 

180 

7 

2 

93 

2 

150 

8 

2 

93 

3 

150 

9 

2 

93 

4 

150 

10 

2 

93 

2 

160 

11 

2 

93 

2 

180 


parts by weight 
parts by weight 


Tensile 

Strength Modulus 

Elonga- 
tion -to 
break 

MPa 

ksi 

GPa 

Msi 

% 

82.1 

11.9 

3.02 

0.44 

5.56 

79.3 

11.5 

2.92 

0.42 

5.22 

82.1 

11.9 

2.90 

0.42 

5.77 

84.8 

12.3 

2.97 

0.43 

6.78 

79.3 

11.5 

2.87 

0.42 

5.36 

77.2 

11.2 

2.88 

0.42 

3.20 

92.4 

13.4 

3.15 

0.46 

7.99 

89.0 

12.9 

3.01 

0.44 

5.89 

91.7 

13.3 

2.93 

0.43 

7.60 

80.0 

11.6 

2.93 

0.43 

5.75 

80.0 

11.6 

2.88 

0.42 

5.55 


*Celanese Plastics and Specialties Co. 
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TABLE XIII. RESIN CASTING PROPERTIES OF SYSTEM "F"* 


RT 


121°C 


Elonga- 

Tensile tioh-to Flexural 


Strength 

Modulus 

break 

Strength 

Modulus 

MPa ksi 

GPa Msi 

% 

MPa ksi 

GPa Msi 


92.4 

13.4 

3.15 

0.456 

7.99 

102.8 

14.9 

2.90 

0.420 

39.3 

5.7 

1.93 

0.280 

9.71 

72.4 

10.5 

2.01 

0.292 


3 Day H-O boil, 97.2 14.1 2.60 0.377 

tested at RT 


Transition Temperature Determined by TMA : 

Initial T g = 176 °C 
After 3 day H 2 0 Boil = 142 °C 
After 7 day H^O Boil = 143 °C 


Heat Deflection Temperature = 155 °C 
C@ 1.82 MPa (264 psi)] 


Notched Izod = 58.7 J/m of notch 

(1.10 ft-lbs/in of notch) 

♦Composition: Epi-Res 508, 100 parts by weight 

Epi-Cure 841, 22.5 parts by weight 

Cure Schedule: 2 hours @ 93 °C 

2 hours § 150°C 
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I 


3. C-6K/ETP 10 -An epoxy -thermoplastic blend (in-house devel- 

oped material) with 10 parts thermoplastic in 
100 parts epoxy. 

4. C-6K/ETP 50 - Same as item 3 except this consisted of 50 

parts thermoplastic in 100 parts epoxy. 

Silicone rubber (a prepolymer) was dissolved in n-hexane, 
then the curing agent was added to the solution before sizing. 
The silicone rubber was mixed at a ratio of 10 parts silicone 
rubber prepolymer to 1 part curing agent. 


The FRE 25 sizing solution was prepared by dissolving the 
organophosphazene prepolymer in acetone. Then, Epon 828 was 
added in the ratio of 25 parts elastomer to 100 parts epoxy. 

Just before sizing, 1% by weight of the elastomer of Esperox 10 
(t-butyl perbenzoate) was added to the solution. This organoper- 
oxide was used as a curing agent for the organophosphazine elas- 
tomer . 

ETP 50 and ETP 10 were dissolved in acetone and the concen- 
tration of all the sizing solutions were adjusted together with 
the sizing speed to allow approximately 1% by weight pickup by 
the fibers . 

Silicone rubber was selected to demonstrate the effect of a 
very soft interphase together with a weak interfacial joint 
strength. FRE 25 was used to illustrate the effect of a second, 
soft and dispersed phase at the interphase. Organophosphazene 
was selected because of its high temperature performance combined 
with a relatively high strength for an elastomeric material ; it 
is also tough and is compatible with epoxy resins. The purpose 
of the polyblend material is to maintain the strength at the inter- 
face while the elastomeric material is damping the impact and/or 
blunting crack tips to enhance toughness . 

The same principles were applied in the selection of the thermo- 
plastic-epoxy blend, ETP 10 and ETP 50. The only difference is 
that the thermoplastic epoxy blends provide higher modulus as well 
as strength. The difference in concentrations is to provide step- 
wise increases in modulus and strength at the interphase. 

Although the properties of the finish variants were not charac- 
terized due to many unsuccessful attempts to make castings from 
these materials, it is estimated that the moduli of the interphases 
should be in the following order: 

ETP 10 > ETP 50 > FRE 25 > SR 
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EXPERIMENTAL RESULTS AND DISCUSSION 
Wetting Force Measurement 

The Wilhelmy micro-balance wetting force technique was used 
to characterize the surface of carbon fibers. A micro -computer 
was used to control the wetting experiment as well as to carry out 
the data processing. The details of the apparatus construction 
are discussed in Appendix B and the program developed by the authors 
for an HP-85A micro-computer is attached in Appendix C. 

To characterize the surface energetics of the carbon fiber, 
double distilled water was used initially. A very careful cleaning 
procedure for the glassware which contains the wetting liquid was 
required so that the surface tension of the contained water approaches 
72.4 x 10“ 3 N/m. This vigorous cleaning process is detailed in 
Appendix D. By employing this cleaning process, a consistent 72.3 x 
10“ 3 N/m ± 1 x 10“ 4 N/m surface tension for water was measured. Con- 
tact angles for unsized fibers were measured successfully with water 
as the wetting liquid. However, when the fibers contained finish, 
misleading wetting forces were obtained. This is primarily due to 
the finish slowly dissolving into and contaminating the water. As 
a consequence, a significant drop in surface tension of thfe water 
took place. This caused the wetting force to drop during the ex- 
periment. Some of these results are listed in Table XIV. The large 
scatter of data within each region with little or no hysteresis 
effects indicates that the surface tension of the wetting liquid 
was changing during the experiment. The surface tension of the 
water dropped as much as 5 x 10“ 3 N/m after experiments were conducted. 

Under this condition, water, as well as other commonly used 
liquids with known surface tension are not applicable for this study 
because of the contamination/dilution process. Nevertheless, in 
order to understand the wetting behavior and to prescribe the optimum 
prepregging condition, a major component of the matrix resin should 
be used as the wetting liquid. Although contamination by finish 
from the sized fiber still continues, this is representative of the 
prepregging process. A DGEBA epoxy resin, Epi-Res 508, was selected 
to be the wetting liquid. Epi-Res 508 was chosen because of its 
relatively low viscosity and it is the major component of the 
matrix used in this study. A viscosity -temperature profile of Epi- 
Res 508 was determined in a Brookfield viscometer and is shown in 
Figure 37. An insulated clean beaker was used to contain and to 
maintain the resin at approximately 70° C. At this temperature, a 
reasonable low viscosity, 0.16 Pa (160 cps) was obtained and no sig- 
nificant dynamic wetting phenomenon was observed. The wetting force 
data obtained for a series of fibers and finishes are shown in 
Table XV. 

Comparing the data in Table XV, the effect of finish variants 
on wetting behavior was demonstrated. The receding angles of all 
the samples are all similarly high suggesting that, once the fibers 
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TABLE XIV. WETTING FORCE MEASUREMENTS OF CARBON FIBERS 
WITH WATER AT 23 °C 


Advancing, yg Equilibrium, yg Receding, yg 

Sample ( contact angle, 9 ) ( contact angle, 6 ) ( contact angle, 9 r 


C-6KU* 1 * 

82.3 

± 

5.5 

89.4 

± 

7.2 

97.5 

+ 

9.7 


(60° 

± 

2.3°) 

(57° 

± 

3°) 

(53° 

± 

4.3°) 

as-i 1 (2) 3 

81.7 

+ 

7.5 

91.3 

± 

9.8 

108.0 

± 

8.0 


(61° 

± 

3°) 

(57° 

+ 

4°) 

(50° 

+ 

3.3°) 

C-6KE ^ 

44.2 

± 

14.6 

41.8 

± 

24.0 

63.5 

± 

10.0 

T-6300 (4) 

90.9 

± 

18.3 

91.3 

± 

23.3 

100.4 

± 

21.4 


(1) Celion® 6000 unsized carbon fiber. 

(2) Hercules AS-1 unsized carbon fiber. 

(3) Celion® 6000 standard epoxy compatible sized fiber, 
commercial product. 

(4) Thornel T-6300 standard commercial sized fiber. 



Pa* s 

(xlO 3 cps) 



Temperature (°< 


FIGURE 37. VISCOSITY-TEMPERATURE PRO! 

EPOXY RESIN MEASURED BY A 
WITH SPINDLES NO. 1 AND 2 



TABLE XV. 


WETTING FORCE MEASUREMENTS OF CARBON FIBERS 
WITH EPI-RES 508 RESIN AT 70 °C 



T-6300 (1) 

72.8 

± 

27.6 

AS-1 (2) 

57.2 

± 

10.6 

C-6KV (3) 

67.0 

± 

11.3 

C-6KU^ 4 * 

75.6 

± 

19.4 

C-6KE (5) 

55.8 

± 

19.5 

C-6K/ETP 50 (6) 

107.4 

± 

28.4 

C-6K/ETP 10 (7) 

86.0 

± 

40.0 

C-6K/FRE 25 (8) 

71.7 

± 

10.7 

C-6K/SR (9) 

-1.6 

± 

18.2 


95.6 ± 18.5 103.8 ± 20.2 

74.4 ± 5.4 87.2 ± 7.8 

80.2 ± 9.9 92.2 ± 14.9 

83.2 ± 19.5 88.6 ± 20.9 

91.4 ± 13.3 102.4 ± 12.9 

127.2 ± 29.4 135.0 ± 29.1 

103.0 ± 39.0 109.0 ± 43.0 

82.0 ± 12.6 90.2 ± 14.9 

50.6 ± 4.1 99.8 ± 48.2 


Note: All carbon fibers were surface treated except C-6KV. 


(1) Thornel standard sized 6000 filament per tow carbon fiber. 

(2) Hercules unsized carbon fiber. 

(3) Celion® 6000 virgin fibers — no surface treatment and unsized. 

(4) Celion® 6000 unsized carbon fiber. 

(5) Celion® 6000 standard epoxy compatible size, commercial product. 

(6) Celion® 6000 sized with finish variant ETP 50. 

(7) Celion® 6000 sized with finish variant ETP 10. 

(8) Celion® 6000 sized with finish variant FRE 25. 

(9) Celion® 6000 sized with finish variant SR. 
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are "wetted," there is little or no rearrangement of the resin 
film on the surface of the fiber. However, in the advancing 
measurements, large variations are observed. This suggests 
that some finish variants have better compatibility with the 
matrix resin and can be wetted easier. For example, silicone 
rubber coated Celion® 6000 showed the lowest wetting force 
(negative) which indicates that, in order to wet this surface, 
additional external forces are needed to forcibly wet it. This 
will most likely also result in poor adhesion. 

On the other hand, the finish variant ETP50 showed the highest 
advancing wetting force which suggests that the resin is highly 
compatible with this surface and it will probably translate to 
good interfacial adhesion. Furthermore, with a high advancing 
wetting force fiber, better quality prepreg is likely, due to 
the minimization of trapped air at the interface by rapid advance- 
ment of the wetting liquid/resin. In other words, improved lam- 
inate properties are anticipated. 

Single Filament Adhesion Test 


As previously discussed, the single filament critical length 
determination was identified to be one of the better methods for 
characterization of carbon fiber/epoxy interfacial bond strength. 
Recall the load transfer mechanism, wherein the critical length, 
l is determined by equation (41) : 


l 

c 



(41) 


where 

a is the tensile strength of the fiber at critical 
length 

d f is the diameter of the fiber 

is the interfacial shear strength 


By determining Cf and measuring £ c and d f , the interfacial shear 
strength, Ty can be obtained as already stated by rearranging 
equation (41) to the form shown for equation (30): 


T 

y 



c 


(30) 
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The critical length and the diameter of the fiber can be measured 
optically by embedding a single filament in an epoxy matrix and 
stressing the specimen to a strain value such that no additional 
fracture can be obtained with additional stress. However, to ob- 
tain the tensile strength of carbon fiber at the critical length 
is itself an involved task. Carbon fiber strength is sensitive 
to gauge length and the critical length of carbon fiber is too small 
to experimentally test in a reliable manner. Nevertheless, it was 
shown (1) that a linear relationship can be obtained by plotting 
single filament tensile strength versus the logarithm of gauge 
lengths. The test method is described in Appendix E. 

Gauge lengths of 7.62 cm (3-inch), 2.54 cm (1-inch) , 1.27 cm 
(%-inch) , 0.635 cm (^-inch) and 0.254 cm (1/10-inch) were tested 
for each fiber type and finish variant and best fit straight lines 
were obtained from these five data points for each variant. These 
straight lines were then extrapolated to the critical length and 
the tensile strengths were obtained. The results of these gauge 
length studies are presented in Appendix F. 

To perform the adhesion test, a single filament of carbon fiber 
was aligned in the center of an ASTM 6.35 cm (2% inch) dogbone epoxy 
casting. The dimensions of the dogbone specimens are 0.32 cm 
(1/8 inch) wide x 0.16 cm (1/16 inch) deep with a 2.54 cm (1 inch) 
long gauge section. Fabrication of the specimen and alignment of 
the fiber was accomplished with the aid of an RTV-630* silicone 
rubber mold. A silicone rubber mold was used because of the ease 
in making a flexible mold cavity by casting the rubber around stand- 
ard ASTM injection-molded dogbone specimens. Moreover, the use of 
silicone rubber molds also prevents fiber surface contamination be- 
cause they do not require the application of mold release agents. 

Sprue spots were cut precisely in the center of each end for 
each dogbone to a depth of 0.08 cm (1/32 inch) to align and secure 
each fiber to be embedded. Once the fiber was aligned, a drop of 
rubber cement was used to hold it in place. During the fiber selec- 
tion, mounting and the application of rubber cement, the section of 
the fibers placed in the mold cavity were not handled nor contamin- 
ated. 


After the rubber cement was set, system "F" resin was intro- 
duced with an eye dropper into the cavity. Extreme care was exer- 
cised to prevent damaging and/or mis-aligning the fiber. Then the 
whole was degassed in a vacuum oven and cured according to the 
schedule previously described. 

Since system "F" is an amber transparent solid, visual and 
microscopic evaluation of the fiber and the critical length can be 
performed with ease. Furthermore, system "F" also exhibits photo- 
elastic properties which can be used to evaluate the stress patterns 
at the interphase . 

*General Electric Co., Silicon Rubber Dept., Waterford, NY. 
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To apply this test to the above mentioned specimen, a special 
tensile fixture was designed and mounted onto a stage micrometer 
for a Leitz ortholux optical microscope. Long working distance ob- 
jectives and polarizing condensers were used. The details of the 
design for the tensile fixture are presented in Appendix G. 

To carry out this study, the specimen was initially examined 
for fiber alignment, breakage and residual stress. Due to the low 
cure and post cure temperatures employed with this system, together 
with a slow cool down rate, no apparent residual stresses were ob- 
served in the samples. The specimen was then stressed to different 
strain values at approximately 0.5% intervals and at each strain 
level, the number of the broken segments was counted. When this 
number remains constant with increasing strain, the critical length 
of the fiber is reached. This test was repeated on at least ten 
specimens for each fiber type/finish variant and the cumulative 
fiber lengths were used to determine the critical length. 

Once the segment lengths of the fibers were obtained, a 
statistical approach to determine the critical length of the fiber 
was used. According to Drzal et al , (2) these data sets fit the 
two parameter Weibull distribution auite well. 

F(X) = 1 - exp for X>o (42) 

where 

X is the aspect ratio, £ c /d^ 

a is the shape parameter 

3 is the scale parameter 

The maximum likelihood methods were employed to estimate the 
parameters of the Weibull distribution. These parameters are 
solutions to 



94 



A computer program was developed to evaluate these two 
parameters, a, f5, numerically, these were then used to calculate 
the shear strength and the listing is attached in Appendix H. 

The mean and the variance of the interfacial shear strength were 
then obtained from 



(45) 




(46) 


where 


T is the 


Gamma function 


Another program was also developed to plot the cumulative 
function and the population density function (histogram) of the 
critical aspect ratio for each fiber type/finish variant. This 
program is also listed in Appendix I. 

The results of the interfacial bond strength determinations 
for the different fiber types/finish variants are listed in 
Table XVI and the histograms of the critical aspect ratios are 
shown in Figures 38 through 46. The results clearly indicate 
that the surface treatment improved the interfacial adhesion 
significantly and that coating the carbon fibers with silicone 
rubber inhibits inter facial bonding. The remainder of the sized 
and unsized fibers together with the other finish variants did 
not show meaningful difference. 

An underlying reason for the lack of additional distinctions 
is the large variation in the shear strengths and the 1:1 rela- 
tionship between extrapolated tensile strength and shear strength. 
This effect is seen in the AS-1 data shown in Appendix F whereby 
an abnormally high extrapolated tensile strength translates 
directly to a proportionally high shear strength. 

In conclusion, this test can be used to identify strong and 
weak interfacial bond strengths successfully. However, due to 
the nature of this test, it cannot reliably differentiate small 
differences . 
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TABLE XVI. SINGLE FILAMENT INTERFACIAL SHEAR STRENGTH 



ff f 

,i c 



T 


Sample 

MPa 

ksi 

a 

6 

MPa 

ksi 

T-6300 (1) (Sized) 

4414 

640 

3.657 

46.368 

60 ± 29 

8.682 + 4.211 

AS-1 ^ (Unsized) 

5931 

860 

3.294 

41.736 

92.6 + 53.2 

13.43 + 7.707 

C-6KV (2) 

2772 

402 

3.255 

71.385 

25.4 + 14.9 

3.687 + 2.159 

C-6KU 

4828 

700 

3.413 

45.282 

68.6 + 37.1 

9.948 + 5.382 

C-6KE ^ 

5379 

780 

3.094 

47.092 

76.3 +48.9 

11.07 + 7.09 

C-6K/ETP 50 

5400 

783 

3.752 

47.394 

71.1 + 33.2 

10.31 + 4.81 

C-6K/ETP 10 

5655 

820 

3.857 

51.726 

67.7 + 30.3 

9.81 ± 4.39 

C-6K/FRE 25 

5172 

750 

3.742 

46.899 

68.9 + 32.3 

9.99 ± 4.68 

C-6K/SR 

4276 

620 

3.001 

185.952 

15.57 ± 10.56 

2.257 + 1.531 


(1) Commercially available carbon fibers 

(2) Virgin carbon fibers with no surface treatment and unsized 
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FIGURE 39. HISTOGRAPH OF THE CRITICAL ASPECT RATIO OF 
AS-1 UNSIZED CARBON FIBER 
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40. HISTOGRAM OF THE CRITICAL ASPECT RATIO OF CELION® 
6000 WITH NO SURFACE TREATMENT AND UNSIZED 
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41. HISTOGRAM OF THE CRITICAL ASPECT RATIO OF CELION® 
6000 WITH SURFACE TREATMENT BUT UNSIZED 
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FIGURE 42. HISTOGRAM OF THE CRITICAL ASPECT RATIO OF CELION® 
6000 STANDARD EPOXY COMPATIBLE SIZE 
(COMMERCIAL PRODUCT) 
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FIGURE 43. 


HISTOGRAM OF THE CRITICAL ASPECT RATIO OF CELION® 
6000 WITH ETP-50 SIZE 


99 



•ET 1 

HE 1 BULL DISTRIBUTION 
P 0 PULATI 0 N D E N S I T Y E T P 1 0 

0 4 0 0 r 


. 0300 

S 0 2 0 0 


W 0 1 0 0 

ct 

u 

0 0000 


ASPECT RA-TIO 


FIGURE 44. HISTOGRAM OF THE CRITICAL ASPECT RATIO OF 
CEL I ON® 6000 WITH ETP-10 SIZE 
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FIGURE 45. HISTOGRAM OF THE CRITICAL ASPECT OF 
CELION® 6000 WITH FRE-25 SIZE 
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FIGURE 46. HISTOGRAM OF THE CRITICAL ASPECT RATIO OF 
CELION® 6000 WITH SILICONE RUBBER SIZE 



Although this test cannot differentiate small differences 
numerically in interfacial bond strength, it does provide some 
information on fracture behavior of different fiber types/finish 
variants. Three distinctly different types of fracture behaviors 
were observed. The optical photomicrographs of these fracture 
behavior patterns are shown in Figures 47 through 49 and a 
schematic diagram of these three fracture modes is presented in 
Figure 50. 

From this information, it is clear that the commercial epoxy- 
based sized fibers provide a higher interfacial shear strength and 
that brittle matrix cracks were initiated at the fiber fracture 
ends. However, ETP-50 sized Celion® 6000 produced a ductile, 
energy absorbing type of crack at the interphase. Further in- 
crease in stress on the latter fibers created a sharp, brittle 
crack extending into the bulk matrix material. These observa- 
tions suggest that it may be possible to achieve both high shear 
strength and increased fracture toughness simultaneously. 
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CROSS POLAR TRANSMISSION PHOTOMICROGRAPH 100X 



TRANSMISSION PHOTOMICROGRAPH 3 2 OX 


FIGURE 47. FRACTURE MODE AT FIBER END OF LOW INTERFACIAL 
BOND STRENGTH (SILICONE RUBBER SIZED AND 
UNSIZED CARBON FIBERS) 




TRANSMISSION PHOTOMICROGRAPH 3 2 OX 



CROSS POLAR PHOTOMICROGRAPH 20 OX 


FIGURE 48. FRACTURE MODE AT FIBER END OF HIGH INTERFACIAL BOND 
STRENGTH WITH STANDARD EPOXY SIZED CARBON FIBER. 
BRITTLE FAILURE INTO MATRIX AT FIBER ENDS. 
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CROSS POLAR PHOTOMICROGRAPH 320X 



TRANSMISSION PHOTOMICROGRAPH 320X 


FIGURE 49. FRACTURE MODE AT FIBER END OF LOWER MODULUS FINISH 
(ETP 50) CARBON FIBER. A DUCTILE FRACTURE MODE AT 
INTERPHASE INDICATES A TOUGH FRACTURE BEHAVIOR. 

THE INNERLAYER IS IDENTIFIED BY THE SHARP BRITTLE 
CRACK EXTENDED INTO THE BULK MATRIX MATERIAL. 
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a) Typical unsized fiber or weak interfacial 
bonding - no matrix crack. 




b) Celion®6000 standard epoxy sized fiber 
brittle cracks initiated at fiber ends. 




c) Celion®6000 sized with ETP-50 

different crack geometry immediately adjacent 
to fiber. Brittle matrix crack developed 
outside the interphase region. 


FIGURE 50. A SCHEMATIC DIAGRAM OF THE THREE TYPES OF FAILURE 
MECHANISMS AT FIBER ENDS WITH DIFFERENT FINISHES 
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Effect of Sizings on Composite Properties 


The finish variants were formulated and applied to Celion® 6000 
unsized carbon fiber. A three-inch wide prepreg tape with a thick- 
ness of 0.127 mm (5 mils) was prepared from each variant. Laminates 
were then fabricated using a match-metal-die molding technique. 
0°-tensile properties together with 0°-flexural, 90°-flexural and 
short beam shear (SBS) at both room temperature and 93.3°C (200°F) 
were measured. An ETI-Dynatup drop dart instrumented impact tester 
(Appendix J) was used to measure the impact toughness of additional 
laminates. An approximately quasi-isotropic , [+45/90/-45/0/ 45/90] s 
lay-up was used for the impact toughness measurements . To character- 
ize the toughess of these laminates, the maximum penetration force, 

Pp recorded on the tester is reported along with Pj , the force assoc- 
iated with initial damage. These composite properties together with 
the wetting force measurements and the single filament adhesion strength 
are presented in Table XVII. 

To further evaluate this toughness enhancement obtained with 
resin system "F," a commercial resin system 5208 is used. Finish 
variant ETP-50 and a standard commercial size (serving as a control) 
were applied to Celion® 6000 unsized fiber and their properties are 
presented in Table XVIII. 

Comparing the wetting force measurements with both the single 
filament adhesion strength and the composite shear strength, it can 
be seen that good wetting is required for good adhesion. However, 
better wetting as measured herein does not necessarily provide stronger 
interfacial bonding. This also suggests that once the fibers are 
"wetted," regardless of how the wetting is achieved, either forcibly 
wet or due to differences in surface energetics, then some adhesion 
can be achieved. This is identified in the silicone rubber sized car- 
bon fiber composites where structural integrity is maintained in spite 
of chemical incompatibility. However, stronger interfacial bond 
strengths are clearly obtained with the other finish variants and the 
control. This probably arises from more favorable chemical bonding 
at the interface. 

Although a 1:1 correlation between the single filament adhesion 
strength and the composite shear strength was not demonstrated, the 
single filament test is still very useful. This test method provides 
valuable insight into adhesion fundamentals, the underlying composite 
shear strength and the failure mechanisms involved in corresponding 
composites . 

It is noteworthy that three of the finish variants selected in 
this study, ETP-50, ETP-10 and FRE-25 all show significant increases 
in impact toughness with no corresponding decreases in other mechani- 
cal properties either at room or elevated temperature. This impact 
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tabu xvii. celio# 6009/ggnw *r unnA-ra. piqronEB 


HBTDC KKX H'BUCT 5 


iawi 

Aoywcsc 

<■91 

wgmjBRim 

AKZDDC 

T-4300 (siaed) 

72.0*27.6 

95.6*10.5 

103.6*20.2 

AS-1 (unalsed) 

57.2*10.6 

74.4*5.4 

07.2*7.0 

0-6 rv 

67.0*11.3 

60.2*9.9 

92.2*14.9 

C-6W 

75.6*19.4 

63.2*19.5 

60.6*20.9 

C-6EI 

55.6*19.5 

91.4113.3 

102.4*12.9 

C-6K/E1W0 

107.4*26.4 

127.2*29.4 

135.0*29.1 

C-6VCTP10 

W. 0*40.0 

103.0*39.0 

109.0*43.0 

.C-6K/FBE-25 

71.7*10.7 

62.0*12.6 

90.2*14.9 

C-4H/B* 

-1.6*11.2 

50.6*4.1 

99.6*40.2 


BDGiz nuwNH tbcilkS 


CRITICAL LltCrTH 

RT 


5373*2 

DBTBMDCS OtM 

tk-U 

9TR 

»a 

(*»M 

PCD 

W*a 

i-i) 

3TR 

«>. 

<*■*) 

PCD 

®a 

J-ii 

59.9*29.0 

(0.66*4.21) 

1669 

(242) 

142.8 

(20.7) 

rat 

NA 

,92.6*53.2 
( 13.43*6.71) 

1669 

(242) 

131.0 

(19,0) 

NA 

w 

25.5*13.9 

(3.69*2.02) 

1641 

(267) 

144.2 

(20.9) 

rat 

NA 

66.6*37.1 
(9.95*5. 36) 

2069 

(300) 

153.0 

122.3 

1952 

(283) 

151.7 

(22.0) 

76.3*48.9 

(11.07*7.09) 

1793 

(260) 

141.4 

(20.5) 

1855 

(269) 

144.8 

(21.0) 

71.1*33.2 

1766 

(256) 

130.3 

(10.9) 

1828 

(265) 

146.2 

(21.2) 

67.7±X.3 

2034 

(295) 

141.4 

(20.5) 

1803 

(273) 

133.1 

(19.3) 

60.9*32.3 

1731 

(251) 

145.5 

(21.1) 

1786 

(259) 

133.1 

(19.3) 

15.57*10.56 

(2.257*1.531) 

1621 

(235) 

127.6 

(18.5) 

1434** 

(208) 

135.9 

(19.7) 


0* rtgagMj 90* FLBCURftL 1 


RT 

93. 3" 

RT 


93.3’C 


PCD 

9TR 


PCD 

am 

«P a 

®a 

"Pa 

"Pa 

Q*a 

"P. 


i~il 



i—ii 

<*»t) 

1800 

121.4 

NA 

99.31 

0.83 

rat 

(261) 

(17.6) 


(14.4) 

(1.20) 


1676 

115.9 

NA 

77.93 

7.17 

NA 

( 243) 

(16.8) 

NA 

(11.3) 

(1.04) 


1331* 

127.5 

NA 

» 

PA 

rat 

(193)* 

(10.5) 





1910 

123.4 

PA 

NA 

rat 

rat 

(277) 

(17.9) 

NA 

rat 

NA 

rat 

1848 

123.4 

1069 

81.30 

7.45 

71.03 

(268) 

(179) 

(155) 

(11.0) 

(1.08) 

(10.3) 

2090 

127.6 

931 

84.83 

7.52 

74.46 

(303) 

(18-5) 

(135) 

M2.3) 

(1.09) 

(10.8) 

1586 

122.1 

979 

66.21 

8.83 

71.03 

(230) 

(17.7 

(142)) 

(9.6) 

(1.20) 

(10.3) 

1931 

127.6 

1090 

74.40 

8.21 

63.45 

(280) 

(18.5) 

(150) 

(10.0) 

(1.19) 

(9.2) 

806.9 

130.3 

745 

26.21 

7.59 

20.00 

(117) 

< 18.9) 

(108) 

(3.0) 

(1.10) 

(2.9) 


as 1 





RT 

93.3'C 

Pi * 

Pf « 

«t< J ) 

(ft-lb) 

"Pa 

(fcel) 

93.10 

(13.5) 

"Pa, 

(tail) 

(lb) 

(lb) 

rat 

NA 

PA 


03.45 

(12.1) 

rat 

PA 

PA 


49.66 

(7.2) 

NA 

PA 

PA 


06.21 

(12.5) 

PA 

PA 

PA 


95.06 

(13.9) 

57.24 

(8.3) 

1180 

(267) 

1401 

(333) 

0.45 

(6.23) 

06.90 

(12.6) 

57.24 

(8.3) 

1468 

( 330) 

1926 

(433) 

10.47 

(7.72) 

85.52 

(12.4) 

56.55 

18.2) 

1472 

(331) 

1939 

(436) 

10.67 

(7.07) 

91.03 

(13.2) 

57.24 

(8.3) 

1397 

(314) 

1837 

(413) 

9.34 

(6.89) 

34.40 

(5.0) 

20.2S 

(4.1) 

1886 

(424) 

2220 

(499) 

11.82 

(8.72) 


* Failed lntralMlnar shear. 

** Palled at tab due to lo* river strength. 

1 Date obtained ae tested. 

2 Data noneallxed to 62t vol. fraction. 

^ fomallsed to O.ISoa 10.060 In.) thickness. 

^ Average of 10 ^vebnene. 

5 Average o ( 30 speclsvns. 

Ptotei lafnect data la an average of 6 samples. Each oompoeite properties la m average of 5 samples. 



TABLE XVIII. CELION® 6000/5208 LAMINATE, PROPERTIES 





TEST 

SAMPLE 

PROPERTY 


CONDITION 

TEMPERATURE 

C-6KE 

C-6K/EPT50 

Impact 1 
P F (N) 



RT 

1348 

1744 

e t (j) 




7.46 

9.53 

Tensile 

STR (MPa) 
MOD (GPa) 


RT 

1669 

140.7 

1834 

142.1 


STR (MPa) 
MOD (GPa) 


132 °C 

1738 

154.5 

1814 

153.1 


STR (MPa) 
MOD (GPa) 

WET 

132°C 

1745 

156.6 

1786 

153.1 

0° Flex 1 (2) 3 

STR (MPa) 
MOD (GPa) 


RT 

1972 

118.6 

2041 

121.4 


STR (MPa) 
MOD (GPa) 


132°C 

1731 

119.3 

1986 

119.3 


STR (MPa) 

WET ' 3 ^ 

132°C 

1076 

869 

90°C Flex 

STR (MPa) 
MOD (GPa) 


RT 

87.6 

9.31 

82.8 

9.03 


STR (MPa) 


132°C 

63.5 

56.6 


STR (MPa) 

WET ^ 

132°C 

31.0 

26.9 

SBS 

STR (MPa) 
STR (MPa) 


RT 

132°C 

117.2 

79.3 

111.7 

80.7 


STR (MPa) 

WET 

132°C 

53.8 

49.0 

±45° Tension 

STR (MPa) 
MOD (GPa) 


RT 

161.4 

4.76 

159.3 

4.35 


STR (MPa) 
MOD (GPa) 


132 °C 

144.8 

3.72 

131.7 

3.24 


STR (MPa) 
MOD (GPa) 

wet* 3) 

132°C 

129.7 

2.69 

112.4 

2.28 


(1) Normalized to 1.5 x 10" 3 meter thickness. 

(2) Normalized to 62% vol. fraction. 

(3) Immersed in water at 71.1°C for 14 days. 

Note: All data except impact are average of 5 specimens. 

An average of 6 specimens is used for impact data. 
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toughness improvement is also shown in the 5208 resin system. The 
mechanisms involved are not fully understood. It is postulated that 
the softer innerlayer at the interphase diminishes some of the in- 
ternal residual stresses and thereby eliminates stress concentrations . 
The significant improvements in impact toughness may arise from damp- 
ing effects of shock waves at each fiber interphaise where an energy 
dissipating material is not present. In any case, this study has 
demonstrated that tailoring of the interphase properties can enable 
significant improvement of composite performance. 
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MODEL DEVELOPMENT 


The analytical model discussed under the section, "Stress 
Transfer from Matrix to Fiber, " is valid only if there is no 
interphase region. However, it was shown by Drzal et al(D and 
by the authors that an interphase does exist and its effect on 
composite properties is significant. Therefore, a modification 
of the existing analytical model is essential. 

From the observations and the results of this study, several 
recommendations can be made by the authors to support the develop- 
ment of a new analytical model which includes a third phase - the 
interphase region. 

I 

Experimentally, the authors have shown that incorporation of 
a lower modulus interphase region results in the observation of 
a different fracture mechanism. This energy absorbent inner layer 
also improves the impact toughness of the composite. Further, 
there is some indication that the transverse matrix crack which 
occurs at a fiber fracture is blunted and that the tensile strength 
of a laminate is slightly enhanced. However, it has to be made 
clear that not only the physical properties of the interphase are 
important, but the chemical nature of the interphase also has to 
be considered. This is because some chemical bonding is required 
to provide strong interfacial shear strength to maximize the load 
transfer across the interface. The thickness of this interphase 
may itself be crucial. An annulor layer which is too thin may 
not be effective in improving impact resistance, while too thick 
a layer may lower other laminate properties such as transverse 
tension and compression, etc. 

A schematic diagram of the modified model with the incorporated 
inner layers is shown in Figure 51. A micromechanistic approach 
can be taken to analyze the stresses in each region. Two approaches 
with different assumptions can be taken. 

First, consider the fiber to be very rigid, when under stress 
and the true deformation of the fiber to be negligible compared 
to the matrix and the inner layer. In this case, a ductile inner 
layer of low shear modulus combines with the matrix to undergo 
large deformations in transferring tensile load to the fiber. 

Or second, consider the fiber to be an inclusion in the inner 
layer material, then treat the fiber/inner layer composite with 
its inherent characteristics as an inclusion in the bulk matrix. 

After the stress states are determined, then classical micro- 
mechanistic approach can be used to study the composite properties. 
The resulting analytical model should include a third phase intro- 
duced by either (or both) of the approaches illustrated in Figure 51. 

To incorporate the experimental results reported here, some 
refining experiments are necessary: 
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U is the shear modulus 


k is the bulk modulus 
and subscript 

m = matrix 
I = interphase 
f = fiber 


FIGURE 51. A SCHEMATIC DIAGRAM OF AN INNERLAYERED INCLUSION 




- What is the modulus of the interphase? 

What is the composition gradient (and property gradient) 
across the interphase? 

What is the effect of interphase thickness? 

What is the degree of bonding required at the fiber 
surface? (Is the perfect bonding assumption a valid one?) 

What is the effect of finish on the load transfer to an 
adjacent fiber? (What is the loaded length of the ad- 
jacent fiber in the vicinity of a fracture?) 

The answers to these questions and additional experimental 
observations would be amenable to analytical modelling to predict 
and account for the observed phenomenon. A valid model must be 
consistent with the observations. 
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CONCLUSIONS AND RECOMMENDATIONS 


An extensive literature review on the role of interface and 
interphase in composite properties has been completed and ex- 
periments were designed to clarify some of the uncertainties 
encountered in this survey. 

Wetting angle measurement on fiber surfaces is the most commonly 
used technique to predict adhesion strength due to wetting. However, 
no direct correlation has been shown between wetting angle measure- 
ments and f iber-to-matrix bond strength. In this study, wetting 
behavior has been characterized by the Wilhelmy wetting force 
measurement and interfacial shear strengths have been studied via 
the single filament critical length determination. The results 
indicate that intimate contact is necessary for adhesion, but 
chemical bonding may be responsible for higher interfacial bond 
strength . 

Further, the results also indicate that by tailoring the inter- 
phase, different failure mechanisms can be achieved without signif- 
icant loss in interfacial bond strength. This approach was in- 
vestigated more thoroughly by an in depth evaluation of laminate 
properties. A significant improvement in impact toughness was 
obtained without loss of static mechanical properties. The con- 
clusions of this study can be summarized as the following: 

1. Good wetting properties are essential to provide good bond- 
ability, but better wetting properties do not necessarily 
provide stronger adhesive joint strengths. 

2. Single filament critical length determination is a good 
analytical tool to characterize carbon fiber/epoxy adhesion 
properties and inherent failure mechanisms. 

3. Single filament adhesion strength and composite shear strength 
are only loosely correlated. 

4. Different failure mechanisms are identified in testing of 
single filaments of carbon fiber with different finish variants 
embedded in an epoxy matrix. 

5. A significant 30% increase in impact toughness was measured 
without loss at room, elevated temperatures and hot -wet 
mechanical properties by incorporation of a " tailored - 
interphase" into the composite. 
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6. A modified analytical model is necessary to better understand 
the effect of the interface and interphase on composite per- 
formance . 

7 . A better understanding of the load transfer mechanisms in- 
volved in this more complex 3-component composite environ- 
ment is necessary to prescribe the optimum interphase proper- 
ties . 
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APPENDIX A. CALCULATIONS FOR THE THICKNESS OF A 1.5% BY 
WEIGHT EPOXY SIZE ON CARBON FIBER 


rr 


Density of carbon fiber, = 1.76 g/cc 

Diameter of carbon fiber, r = 7 x 10~4 cm 

Density of epoxy size, P_ = 1.182 g/cc 

s 

Assume unit length, H = 1 cm 

Volume of the fiber, Vf = = 1.539 x 10“6 cc 

Weight of the fiber, Wf = VfPf = 2.709 x 10“6 g 
.•. Weight of the eooxy size, W s = 1.5% Wf 

= 4.064 x 10 -8 g 

Volume of the epoxy size, V s = W s /p s 

= 3.438 x 10“ 8 cc 

But, Vg = ir(R^-r2)£ = 3.438 x 10~8 cc 
.*. R = 7.08 x 10~4 cm 
Thickness of epoxy size, A = R-r = 800 A 
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APPENDIX B 


APPARATUS FOR WETTING FORCE MEASUREMENT 


A Cahn 26 automatic electrobalance was used to measure the 
wetting force of a single filament carbon fiber. This electro- 
balance was interfaced through a BCD interface module to an 
HP-85A microcomputer. This arrangement allowed the data to be 
stored, plotted and manipulated during the experiment. 

In order to assess the advancing, receding and equilibrium wetting 
force, a constant speed moving stage for the wetting liquid 
was required. This stage was made from a modified micromanipulator 
and a high torque, variable speed electric motor. The electric 
motor was controlled automatically by the microcomputer to allow 
automation of this experiment. The program developed for this 
process is listed in Appendix C and a schematic diagram of this 
set-up is presented in Figure B-l. 

A single filament of carbon fiber was selected and attached to 
the end of a thin gauge hang-down wire with rubber cement. This 
fiber was then connected to the balance and the weight of the 
fiber was tared. The wetting liquid contained in a perfectly 
cleaned beaker was raised by the micromanipyilator until the fiber 
end was just immersed in the liquid. Then the program would 
automatically accumulate the data from the microbalance and control 
the stage to obtain the wetting force during the advancing, station- 
ary and receding positions. 
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Closed Chamber to 
Eliminate Air Current and Dust 


FIGURE B-l. APPARATUS ASSEMBLY FOR WETTING FORCE MEASUREMENT 





APPENDIX C 


WETTING FORCE MEASUREMENT PROGRAM 


10 REM >AUT0ST 

15 OPTION BASE 1 
40 ON ERROR GOTO 80 
50 LOADBIN 'REDZER' 

60 OFF ERROR 
70 CHAIN 'WETBAL' 

80 IF ERRN=25 THEN 100 

90 PRINT ■Error number ■ »ERRN» * was encountered at line ■ ? ERRLJ ' » ■ 8 GOT 
0 120 

100 PRINT 'A binary program is already in memory causing an error* Plea 
se * 

105 PRINT 'eliminate this binary (i»e* by turning machine off or execut 
ins * 

110 PRINT ''SCRATCH')* Then you may begin the program again*' 

120 END 
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Introduction 


REDZER is a binary program designed to enhance array manipulation capa- 
bilities of the HP-85. An array can be reorganized (redimensioned) or 
initialized to zero with use of one of the two statements offered. This 
binary program can be loaded into the HP-85 by the following command: 

LOADBIN "REDZER" 

For a further discussion of binary programs, see page 193 of the HP-85 
Owner's Manual and Programing Guide. 

ARRAY (RE ORGANIZATION 

You can reorganize an array into a more useful configuration by redim- 
sioning that array. This changes the working size of the array; subse- 
quent statements affect only the elements included in the new working 
size. However, elements not included in the new working size are still 
associated with the array. The values of these elements are not changed, 
and they can be accessed if the array is again redimensioned. 


REDIM array (redim subscripts) [.array (redim subscripts)...] 


The redimensioning subscripts are numeric expressions, variables, or con- 
stants that specify a new upper bound for each dimension. The number 
of subscripts must be the same as the number specified in the original 
DIM, REAL, SHORT, or INTEGER statements. Furthermore, the total number 
of elements in the new working size cannot exceed the number originally 
dimensioned. 
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Examples: 

DIM A(4) ,B(2 ,4) ,C(1 ,9) 

REDIM A(3) 

REDIM B(3,2) 

REDIM A(4),B(2,4) 

X=2 @ REDIM C(2*X-1,10/X-1) 


OPTION BASE 0 assumed. 

Redimensions working size from 
five to four elements. 

Redimensions B from 3x5 matrix 
(15 elements) into 4x3 matrix 
(12 elements). 

Redimensions A and B back to 
original sizes. 

Redimensions C from 2 x 10 matrix 
into 4x5 matrix. 


When a matrix is redimensioned, the values of its elements are reassigned 
to different positions within the matrix. Values of matrix elements are 
stored in order from left to right along each row, from the first row 
to the last. The redimensioning takes the elements out of the matrix 
in that order, and reassigns them in accordance with the new working 
size of the matrix. 

The following example shows how values of matrix elements are reassigned 
when a matrix originally declared to be 3 x 3 is redimensioned into a 
2x2 matrix. The values of the original matrix elements are integers 
that indicate the order in which the elements are stored before the re- 
dime nsioni ng. 

Example: 

10 OPTION BASE 1 
20 DIM A(3,3) 

30 DATA 1, 2,3,4, 5,6,7 ,8, 9 
40 FOR 1-1 TO 3 
50 READ A(I ,1) ,A(I ,2) ,A(I ,3) 

60 NEXT I 

70 K=3 0 GOSUB 130 
80 REDIM A (2 ,2) 


Array A is originally dimensioned as 
3x3. 


Assigns values 1 through 9 to the 
elements of matrix A. 


Displays original 3x3 matrix. 

Redimensions A down to a 2 x 2 
matrix. 
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Ill III 


I II I 


90 K=2 @ GOSUB 130 
100 REDIM A(3,3) 

110 X=3 @ GOSUB 130 
120 END 

130 FOR 1=1 TO X 
140 FOR J=1 TO X 
150DISP A (I ,J) ; 
160 NEXT J 
170 DISP 
180 NEXT I 
190 DISP 
200 RETURN 


Displays redimensioned 2x2 
matrix. 

Redimensions A back up to a 
3x3 matrix. 

Displays redimensioned 3x3 
matrix. 




Subroutine that displays matrix A. 




rmn 

123 Nine elements of original 3x3 

4 5 6 matrix - 

7 8 9 


1 2 

3 4 

12 3 

4 5 6 
7 8 9 


Four elements of redimensioned 
2x2 matrix. Values of elements 
have been reassigned sequentially 
within new working size. 

Nine elements of 3 x 3 matrix with 
original values still assigned. 


Note that redimensioning a matrix does not isolate a submatrix. In other 
words, if you redimension a 3 x 3 matrix into a 2 x 2 matrix, the resulting 
matrix is not the 2x2 submatrix from the upper left corner of the original 
matrix. 
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5 ! WILHELMY BALANCE BATA ACQUISITION .. .PROGRAM WETBAL ♦ , . G . V . NELSON 5/7 
/8 1 . . .VERSION 5.12.81 

8 OPTION BASE 1 

9 DIM D(lfl400)fS2<10)fSl$C60]?Vl*E843 

10 DIM Y(700>rM(3>fS<3)f A*C203 

15 1=0 e 02=1400 

16 DISP USING '7/' i 8 DISP ■ WILHELMY BALANCE' 8 DISP ■ D 

ATA ACQUISITION' 

17 DISP 'Data fro* tape file(Y/N)' 8 INPUT Y* 

18 IF UPC*(Y*C1f13)='Y* THEN 2000 

21 DISP USING '3/' i 8 DISP 'Enter data collection interval (1-3 seconds 
/point ) ' 

22 INPUT S2 

23 DISP 'Enter no. of points to be collected for Advancing Region <300 m 

3 X ) * 

24 INPUT N1 

26 DISP 'Enter no. of points to be collected for Still Region (200 max)' 

27 INPUT N2 

28 CLEAR 8 ON KEY* 1* 'START' GOTO 50 

29 ON KEY* 6 r * UP * GOSUB 6 200 

30 ON KEY* 7 f ' DOWN* GOSUB 6250 

31 ON KEY* 8 . * STOP ' GOSUB 6300 

32 OFF KEY* 3 8 OFF KEY* 4 8 OFF KEY* 2 

37 ! 

38 GCLEAR 8 CLEAR 8 KEY LABEL 

39 GOSUB 6000 ! INITIALIZE INTERFACE 

40 GOTO 40 

50 DISP 'Enter Sample Description (20 Chars max* no commas)* 

51 INPUT A$ 

52 PRINT USING '7/' J 8 PRINT ' WILHELMY BALANCE' 8 PRINT ' 

DATA ACQUISITION' e PRINT 

53 PRINT A$ 

54 PRINT 8 PRINT 'DATA INTERVAL=* . S2 . * SECONDS/POINT' 

55 PRINT 'ADVANCING REGION* ' . N1 . ' POINTS' 

56 PRINT 'STILL REGION* ' » N2 » ' POINTS' 

57 PRINT 'RECEEDING REG ION= ■ t N1 i ■ POINTS' 

58 PRINT 8 PRINT 

80 Y 1 =0 8 Y2=0 8 M=0 

85 GOSUB 7000 8 D1=D/100 8 D2=D1 ! START TIMER 

86 LOCATE 0.130.0.100 

87 FRAME 

90 SCALE 0. 2*N1+N2,0,D1 
100 MOVE NlrO 8 DRAW Nl.Dl 

102 MOVE N1+N2.0 8 DRAW N1+N2.D1 

103 LORG 6 8 MOVE N1/2.D1 e LABEL 'ADVANCING' 

104 MOVE N1+N2/2.D1 8 LABEL 'STILL' 

105 MOVE N1+N2+N1/2.D1 8 LABEL 'RECEDING' 

120 MOVE 0.0 

135 GOSUB 6230 ! MOTOR ON-FWD 

140 1=0 8 J=I 

150 IF I=J THEN 150 

155 GOSUB 5030 

160 Y< I ) . Y=D 

170 J=I 

180 GOSUB 770 

182 PEN 1 

185 MOVE I.Y(I) 
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190 PLOT I r Y ( I ) 

200 IF I=N1 THEN GOTO 220 
210 GOTO 150 
220 ! 

225 GOSUB 6330 ! MOTOR OFF 
230 PRINT 

250 OFF TIMER# 1 8 PRINT * ADVANCING* * 

254 PRINT USING * 10A » 6D ♦ 3D r / » 10A » 6D . 3D ■ * ■ MEAN = 'JM»' STDEV = 'JSQ 
R < SI ) 

260 PRINT 

270 M ( 1 ) =M 8 M=0 8 Y1=0 8 Y2=0 

280 ! TURN ON TIMER FOR STILL REGION 

290 GOSUB 7000 ! START TIMER 

310 1=0 e J=I 

320 IF I=J THEN 320 

325 GOSUB 5030 

330 Y(N1+I ) f Y=D 8 J=I 

340 GOSUB 770 

342 MOVE N1+I»Y<N1+I) 

344 DRAW Nl + I » Y < N1 + 1 ) 

346 IF I =N2 THEN 390 
350 GOTO 320 

390 OFF TIMER* 1 8 PRINT 'STILL * * 

400 PRINT USING ■ 1 OA » 6D . 3D » / 1 10A » 6D * 3D * t ■ MEAN = 'iHf 1 STDEV = '»SQ 
R ( SI ) 

410 PRINT 

420 M ( 2 ) =M 8 M=0 8 Y1=0 8 Y2=0 

430 ! TURNON TIMER FOR RECEDING REGION 

450 GOSUB 7000 ! START TIMER 

455 GOSUB 6280 ! MOTOR REV 

460 1=0 8 J=I 

470 IF I=J THEN 470 

475 GOSUB 5030 

480 Y ( N1 + N2+ 1 ) f Y=D 8 J=I 

490 GOSUB 770 

492 MOVE Nl+N2+IrY( N1+N2+ I ) 

494 DRAW N1+N2+I »Y(N1+N2+I ) 

496 IF I=N 1 THEN 510 

500 GOTO 470 

510 ! STOP ROUTINE 

520 GOSUB 6300 ! STOP MOTOR 

530 OFF TIMER# 1 8 PRINT 'RECEDING REGION? ' 

540 PRINT USING ' 10A » 6D . 3D » / » 10A r 6D . 3D ' » • MEAN = 'JMJ* STDEV = 'JSQ 
R ( SI ) 

550 PRINT e COPY 

554 ON KEY* 5 1 'STORE ' GOTO 1000 

555 ON KEY* 3»'LIST' GOTO 900 

556 ON KEY* 4»'PL0T' GOTO 560 

557 GOTO 37 

558 ON KEY* 2 » ' READ* GOTO 1900 

559 GOTO 37 

560 \ REPLOT GRAPH TO SCALE 

561 GCLEAR 

562 DISP 'WHICH REGION DO YOU WANT TO PLOT ( A=ADVANCING r S=STILLJ R=RECE 
DING) ' 

563 INPUT Al*8 N3=0 

564 IF UPC*(A1*C1»13)='A' THEN N3=l 8 N4=N1 8 M$=' ADVANCING' 

565 IF UPC*(A1*C1»13)='S' THEN N3=N1+1 8 N4=N1+N2 8 M*='STILL* 
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I 


56A IF UPC*<A1*C1f1:3> = *R* THEN N3=N1+N2+1 0 N4=2*N1+N2 0 M*= * RECEDING * 

567 IF N3=0 THEN 562 

570 PRINT USING '4/* 9 

580 Y8=- 1 « E99 0 Y9=1«E99 

590 FOR I=N3 TO HA 

600 Y8=MAX <Y(I)fY8> 0 Y9=MIN ( Y ( I ) f Y9 ) 

610 NEXT I 

611 CLEAR 0 DISP * MIN Y VALUE='fY9 0 DISP * MAX Y VALUE=*fY8 0 DISP 'DO Y 
OU WANT TO SET SCALE ( Y/N > * f 

612 INPUT A1$0 IF UPC*(A1*C1 rlDs'Y' THEN 613 ELSE 614 

613 DISP • ENTER Ymin^Ymax* 0 INPUT Y9 fY80 GOTO 615 

614 Y9=0 0 Y8=D 1 

615 LOCATE 30f120f20f80 

616 SCALE N3-1 » N4 f Y9 f Y8 

617 FXD 0 1 3 0 LAXES -(N4-N3 + 1 >/10f < Y8-Y9 )./10 f N3-1 fY9f2f2 

618 SETGU 0 MOVE 70f3 0 LORG 4 0 LABEL 'TIME<sec>' 

619 MOVE 8 » 45 0 LDIR 90 0 LABEL *F0RCE(mS>* 0 LDIR 0 0 SETUU 
635 MOVE (N3+N4)/2 fY8+(Y8-Y9)#3/20 0 LORG 4 0 LABEL M* 

640 MOVE N3-1 t 0 
650 FOR I =N3 TO N4 
655 MOVE I f Y < I ) 

660 LORG 5 

665 IF N4-N3>40 THEN PLOT IfY(I> ELSE LABEL *+* 

670 NEXT I 
680 COPY 

690 PRINT USING *7/* 9 

715 GOTO 37 
720 END 

770 ! RUNNING AVERAGE ? STD DEV 

780 M=< ( 1-1 >*M+Y>/I 

790 Y1=Y1+Y 

800 Y2=Y2+Y*Y 

810 S1MY2-Y1*Y1/I)/I 

820 RETURN 

900 ! LIST 

902 CLEAR 

904 DISP 'WHICH REGION DO YOU WANT TO LIST < A=ADVANCING f S=STILLf R=RECE 
DING) ■ 

906 INPUT A1*0 N3=0 

908 IF UPC*<A1*C1f 13)='A* THEN N3=l 0 N4=N1 0 M*=' ADVANCING* 

910 IF UPC* (A1*C1f13)=*S* THEN N3=N1+1 0 N4=N1+N2 0 M*='STILL' 

912 IF UPC*<A1$C1f13)=*R* THEN N3=N1+N2+1 0 N4=2*N1+N2 0 M*=*RECEDING* 

914 IF N3=0 THEN 904 

915 PRINT USING *3/* 9 0 PRINT M* 0 PRINT *TIME<sec) FORCE (m3) * 

920 DISP * ENTER NUMBER OF DATA POINTS PER PRINTED VALUE' 0 INPUT N5 
950 FOR I-N3 TO N4 STEP N5 

952 F9=3 

954 IF D2>20 THEN F9=2 
956 IF D 2 < = 2 THEN F9=4 

958 ON F9 GOTO 959 f 959 f 960 f 965 

959 PRINT USING *DDDZ«DDf4XfDDDZ«DD* 9 I*S2fY(I> 0 GOTO 970 

960 PRINT USING * DDDZ ♦ DD » 4X » DDDZ • DDD * 9 I*S2fY(I> 0 GOTO 970 
965 PRINT USING * DDDZ . DD f 4X f DDDZ ♦ DDDD * 9 I*S2fY(I> 

970 NEXT I 

980 PRINT USING *7/* 9 
990 GOTO 37 

1000 ! STORE ACCORDING TO BSDM FORMAT 
1010 01=2*N1+N2 0 N7=2 0 S7=3 
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1020 V1*=*TIME FORCE ■ 

1030 S1*=*ADV* STILL RECED.* 

1040 S2 ( 1 > =1 0 S2 ( 2 ) =N1 +1 0 S2(3)=N1+N2+1 
1045 REDIM D< N7 »01 ) 

1050 FOR 1=1 TO 01 
1060 D ( 1 r I ) =S2# I 
1070 D( 2 » I > =Y ( I ) 

1080 NEXT I 
1140 CLEAR 

1150 DISP ' * * * * STORE DATA * * * * * 0 PRINT 

1160 DISP 'File names must be <= 6 charac- ters long* (Type to exit) 

a 

1170 DISP "Name of data file = *J0 INPUT F* 

1180 IF LEN ( F$ ) >0 AND LEN(F$><7 THEN 1210 
1185 IF F*= * E * OR NOT LEN ( F$ ) THEN 1530 
1190 GOSUB 1590 
1200 GOTO 1140 

1205 IF F$C1»11=*E* AND LEN(F*>=1 THEN 1530 

1210 DISP ‘Is data medium placed in tape drive(Type 'E' to exit) 1 ?© IN 
PUT N$ 

1220 ON FNA ( N$ > GOTO 1210* 1230* 1210* 1530 
1230 CLEAR 

1240 ON ERROR GOSUB 1710 
1250 X=0 

1260 ASSIGN* 1 TO F* 

1270 OFF ERROR 

1280 IF X=1 THEN 1360 

1290 DISP ‘A file of this name has been found* Do you wish to store* 

1300 DISP 'present data set under this name and destroy old data set 

■*0 INPUT N$ 

1310 ON FNACN*) GOTO 1290*1370*1140*1530 

1330 GOSUB 1590 

1340 CLEAR 0 GOTO 1290 

1360 CREATE F$* (8+02*8) DIV 700+2*700 

1370 ASSIGN* 1 TO F* 

1375 CLEAR 

1390 PRINT* 1*1 * A* f 01 * N7 * VIS * S7 » SI* * S2 ( * > 

1400 READ* 1*2 
1410 PRINT* 1 i D(») 

1420 ASSIGN* 1 TO * 

1430 CLEAR 

1440 PRINTER IS 2 

1450 PRINT ********************************** J •*• JTAB(32) » *** 

1460 PRINT *** » TAB ( 14) f 'STORE* »TAB(32) f * * * 

1470 PRINT *** »TAB(32) »***» ********************************** 

1480 PRINT 

1490 PRINT 'Data and related information is stored in *JF$}*»* 

1500 PRINT 

1510 DISP * Is program medium replaced in device *50 INPUT N$ 

1520 ON FNA ( N* ) GOTO 151 0 > 1 530 f 15 10 t 1530 
1530 CLEAR 
1560 GOTO 37 
1590 BEEP 

1600 DISP 'Improper response — please try again** 

1610 WAIT 1500 
1620 RETURN 

1710 IF ERRN=67 THEN X=1 
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1720 RETURN 
1730 DEF FNA ( At > 

1740 IF LEN<At>=0 THEN FNA=1 ELSE FNA=P0S< ■ YNE" t UPCt < AtC 1 » 1 3 ) ) +1 
1750 FN END 

1900 ! TRAP FOR UNINTENDED LOSS OF DATA 
1905 CLEAR 

1910 DISP •Reading data fro* tape will overwrite current memory* 

1920 DISP ’Continue (E to exit>’»0 INPUT Nt 
1930 ON FNA<Nt> GOTO 1920 f 2000 » 554 f 554 
2000 ! READ DATA FROM TAPE BSDM FORMAT 
2005 CLEAR 

2010 DISP 'Name of data file=*F0 INPUT Ft 

2020 DISP ’Is data medium placed in tape drive(T«pe 'E' to exit>*J0 IN 

PUT Nt 

2030 ON FNA(Nt) GOTO 2020 t 2040 f 2020 »2190 
2040 CLEAR 

2050 ON ERROR GOTO 2080 
2060 ASSIGN# 1 TO Ft 
2070 OFF ERROR 0 GOTO 2090 

2080 DISP ’File ’i Ft?’ not found' 0 GOTO 2010 
2090 READ# IfI i At»01»N7»Vlt»S7fSlt»S2(») 

2100 REDIM D ( N7 f 01 > 

2110 READ* 1 f 2 
2120 READ* 1 } D(f> 

2130 ASSIGN* 1 TO * 

2140 PRINT ’DATA ENTERED FROM TAPE FILE Ft 0 PRINT 
2150 N1 =S2 ( 2 ) - 1 0 N2=S2(3)-S2(2) 

2160 FOR 1=1 TO 01 
2170 Y < I ) =D< 2 » I ) 

2180 NEXT I 

2185 S2=D<1f2)-D<1f1> 

2186 D1 =D ( 2 f01/2 ) 

2190 GOTO 554 
4900 ! 

4910 ! INTERRUPT PROCESSING 
4920 ! 

5000 ENTER 3 i G3»G4 ! READ VALUE 

5001 IF ABS ( Q3 > <60000 THEN 5000 
5010 1=1+1 

5020 RETURN 

5021 ! 

5022 » INPUT DECODING 

5023 ! 

5030 Q1=Q3 

5031 Q2=Q4 

5032 IF Q1*Q3 THEN 5030 

5040 IF INT < Q2/10 > = 1 THEN 5045 
5042 DISP •Overrange* 

5045 Q2=RMD ( Q2 f 10 ) 

5050 IF Q2>0 AND Q2<9 THEN 5100 
5060 ! 

5070 02=1 

5075 DISP ’Exponent Error’ 

5080 GOTO 5190 

5090 ON Q2 GOTO 5140t 5160»5060»5180» 5060 » 5060 » 5060 » 5185 
5100 ON Q2 GOTO 5140 r 5160 t 5060 » 5180 » 5060 » 5060 » 5060 » 5185 
5140 Q2= • 1 
5150 GOTO 5190 
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5160 Q2= ♦ 01 
5170 GOTO 5190 
5180 02= ,001 

5184 GOTO 5190 

5185 02= » 0001 
5190 D=Q1*Q2 
5210 RETURN 
6000 ! 

6010 ! DATA INITIALIZATION 
6020 ! 


6030 

RESET 3 





6040 

CONTROL 

3*1 

• 

9 

0 ! 

NO INTERRUPTS 

6050 

CONTROL 

3*2 

* 

9 

0 ! 

CLEAR OUTPUT BITS 

6060 

CONTROL 

3*3 

A 

9 

5 ! 

5 DIGITS 

6070 

CONTROL 

3*4 

A 

9 

0 ! 

NO EXPONENT 

6080 

CONTROL 

3*5 

9 

2 ! 

2 FUNCTIONS 

6090 

CONTROL 

3*6 

A 

9 

0 ! 

NO DECIMAL 

6100 

CONTROL 

3*7 

* 

f 

0 ! 

NO HANDSHAKE 

6110 

CONTROL 

3*8 

» 

9 

0 


6120 

CONTROL 

3*9 

A 

9 

15 


6130 

CONTROL 

3*10 

A 

9 

0 


6140 

RETURN 





6200 

( 





6210 

! MOTOR 

ON-FWD 


6220 

! 





6230 

CONTROL 

3*2 

A 

9 

1 


6240 

RETURN 





6250 

! 





6260 

! MOTOR 

ON-REV 


6270 

! 





6280 

CONTROL 

3*2 

9 

3 


6290 

RETURN 





6300 

i 





6310 

! MOTOR 

OFF 




6320 

i 





6330 

CONTROL 

3*2 

9 

0 


6340 

RETURN 





7000 

GOSUB 5000 ! 

READ 

START VALUE 


7005 GOSUB 5030 ! DECODE VALUE 
7010 D=49999*Q2 ! SET MAX VALUE 
7020 ON TIMER* 1*1000*S2 GOSUB 5000 
7030 RETURN 
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APPENDIX D. PROCEDURES FOR CLEAN GLASSWARE TO PREVENT 
CONTAMINATION OF WETTING LIQUIDS 


1. Glassware is placed in boiling 1:1 diluted nitric acid for 
five minutes. 

2. Then it is rinsed with double distilled water. 

3. It is next steamed with double distilled water for thirty 
minutes . 

4. It is then vacuum dried at 150 °C for ten hours. 

5. Finally, it is individually wrapped in aluminum foil and 
stored in a dessicator. 
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APPENDIX E 


SINGLE FILAMENT GRAPHITE FIBER TENSILE TEST 


A representative number (generally 20) of single filaments 
are randomly selected from a tow sample. The filaments are center- 
line mounted on special slotted paper tabs as shown in Figure E-l. 
The paper tabs are gripped such that the test specimen is uniaxially 
aligned in the jaws of an Instron. An Instron constant cross-head 
speed of 0.02 inch/min and a chart speed of ten inches/min are 
employed. The procedure of mounting the graphite fiber on the tabs 
is as follows: 

1. The strand bundle is placed loosely on a suitable work surface. 

2. A dental pick is used to gently separate the filaments and a 
suitable test specimen is selected. 

3. A tab of pre-punched slot is used as the mounting medium. 

The single filament is centered along the slot of a mounting 
tab and a small strip of tape is used to anchor one end to 
this tab. 

4. The filament is lightly stretched and the opposite end is 
anchored in the same manner . 

5. A small amount of sealing wax is carefully placed on top of 
the filament at each edge of the slot to secure the filament 
to the tab. 


E-l 




Cardboard Tab 
Wax Bead 


Graphite Fiber 


Cut Tab Here 


Masking Tape 


figure E-l. graphite fiber mounted on a tab 



APPENDIX F. SINGLE FILAMENT TENSILE STRENGTH 
AT VARIOUS GAUGE LENGTHS 


F-l 


MPa 





10 ”« 


icr J 


10" z 


0.1 0.25 0.5 1 


(2.54x10 4 cm) (2.54xl0“ 3 an) (2.54xl0“2an) (0.254an) 

Gauge Length (in) 


(2.54an) 


FIGURE F-2 . AS-1 (STANDARD UNSIZED PRODUCT) LOT #156-3 
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io _J 


io -2 


O.l 0.25 0.5 1 


(2.54xl0 _4 an) 


(2.54xl0 -3 an) (2.54xl0- 2 an) (0.254an) 

Gauge Length (in) 


(2.54cm) 


FIGURE F-3 . C-6KV (UNSIZED. NONTREATED) LOT #7411 






MPa 





MPa 



]_0 - 4 io - 3 10 -2 0.1 0.25 0. 

(2.54xl0- 4 an) (2.54xl0 -3 an) (2.54xl0 -2 an) (0.254cm) 

Gauge Length (in) 


5 1 

(2.54cm) 


FIGURE F-6 . C-6K/ETP-50 LOT #9431 
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10 _J 


10-2 


0.1 0.25 0.5 1 


(2.54xl0~ 4 an) (2.54XlO _ 3an) (2.54x10 2cm) (0.254am) 

Gauge Length (in) 


(2.54am) 


FIGURE P-7. C-6K/ETP-10 LOT #9Y31 





MPa 


Pm 



(2,54xl0-4cin). (2.54xl0-3on)_ 


10-2 0.1 0.25 0 

( 2 . 54xl0 - 2am) C0,254an) 

Gauge Length (in) 


i 

VO 


5 1 
(2,54cm) 


FIGURE F-8. C-6K/F RE-25 LOT #9732 






id p» o> tr cr cr fupi p» 


APPENDIX G. TENSILE FIXTURE TO STUDY SINGLE FILAMENT 
ADHESION PROPERTIES 



b=9 . 5xl0 -3 m 
c=l. 27x10 3 m 
d=2 . 54xl0 -2 m 
e=l . 9xl0 -2 m 
f =6. 4x10 3 m 
g=7 .9xl0" 3 m 
h=l . 3xl0 -2 m 
i=8xl0 _4 m 


FIGURE G-l. TENSILE FIXTURE FOR LEITZ OPTICAL MICROSCOPE 
DRAWN TO SIZE 
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APPENDIX H. 


SINGLE FILAMENT INTERFACIAL SHEAR 
STRENGTH DETERMINATION PROGRAM 


The purpose of this program is to apply a Weibull Distribution 
statistical analysis of the critical aspect ratio of the fibers 
where Weibull parameters as well as the average shear strength and 
the variance of the shear strength are determined. 


H-l 



10 !’*** INTERFACIAL SHEAR STRENGTH PROGRAM -- CRITICAL LENGTH DETERMINA 
TION *** 

20 DIM D ( 500 ) » L ( 500 ) * R ( 350 ) 

30 ON KEY# 1» 'ENTER' GOTO 170 
40 ON KEY* 2? 'STORE* GOTO 600 
50 ON KEY* 3 » * DISPLAY * GOTO 760 
60 ON KEY* 4 » ‘ EDIT * GOTO 940 
70 ON KEY# 5 » * CALC * GOTO 1200 
80 ! ON KEY* 6»‘HIST0* GOTO 
90 ! ON KEY* 7r ‘WEIBULL* GOTO 

100 ! ON KEY# 8 1 ‘ NORMAL • GOTO 
110 ! 

120 CLEAR 8 KEY LABEL 

130 DISP * INTERFACIAL SHEAR STRENGTH DETERMINATION!* 

140 DISP 'CRITICAL LENGTH MEASUREMENT.* 

150 DISP ‘SELECT OPTION* 

160 GOTO 160 

170 ! **** DATA ENTRY **** 

180 CLEAR 8 DISP ‘ENTER DATA FROM KEYBOARD (K) OR TAPE (T)?* 

190 INPUT Y ♦ 

200 IF Y*= * K * THEN 250 
210 IF Y$= * T * THEN 500 
220 GOTO 180 
230 ! 

240 ! 

250 ! *** ENTER DATA FROM KEYBOARD *** 

260 ! 

270 CLEAR 8 DISP 'ENTER SAMPLE NAME?* 

280 INPUT A$ 

290 DISP * MEAN VALUE OF FIBER CRITICAL LENGTH (MICRON) ISJ*»8 INPUT W1 
320 DISP ‘TENSILE STRENGTH (PSI) OF FIBER AT‘?Wlf GAGE LENGTH IS{‘»8 IN 

PUT S 

350 DISP ‘THE TOTAL NUMBER OF SPECIMEN TESTED 4 . *»e INPUT N1 

380 DISP ‘THE TOTAL NUMBER OF BREAKS } ‘ f 8 INPUT N 

410 DISP ‘ENTER FIBER DIAMETER AND CRITICAL LENGTH (MICRON)* 

420 FOR 1=1 TO N 

430 DISP *D( * » 1 1 * ) » L(*»I»*) = ‘ 

450 INPUT D(I>fL(I) 

460 NEXT I 
470 GOTO 120 
480 ! 

490 ! *#** DATA ENTRY **** 

500 ! 

505 DISP * INSERT DATA FILE TAPE! (WAIT FOR THE NEXT COMMAND! * 

506 BEEP 

507 WAIT 3000 

510 DISP ’ENTER FILE NAME 4 * » 8 INPUT A$ 

520 ASSIGN* 1 TO UPC*(A$) 

530 READ* 1 i A$»Wl»SfNl»N*D( > »L( ) 

570 ASSIGN* 1 TO * 

580 GOTO 120 
590 ! 

600 ! *** STORE DATA ON TAPE *** 

610 ! 

620 DISP 'INSERT DATA FILE TAPE! CWAIT FOR THE NEXT COMMAND! 3* 

621 BEEP 

622 WAIT 3000 

630 CLEAR 8 DISP * ENTER FILE NAME * » 8 INPUT A* 
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640 DISP "NEW OR EXISTING FILE <N/E>?’ 

650 INPUT Q$0 IF UPC*(Q*[:i 1 13 > = ’N’ THEN 660 ELSE GOTO 670 
660 CREATE UPC*(A*)n50 
670 ASSIGN# 1 TO UPC#(A*> 

680 PRINT# 1 t A$nW1»SnN1nNnD< ) »L( > 

720 ASSIGN* 1 TO * 

730 DISP ’ DAT A ON • i UPC* ( A* > » • STORED ON TAPE* 0 KEY LABEL 
735 WAIT 3000 
740 GOTO 120 
750 ! 

760 ! *** DISPLAY DATA *** 

770 ! 

780 CLEAR 0 DISP ■ DATA IN MEMORY? <Y OR N>’»0 INPUT T$ 

790 IF T$=’Y* THEN 810 
800 IF T *= ’ N * THEN 120 

805 DISP •ENTER Y OR N’ 

806 GOTO 780 

810 CLEAR 0 PRINT “SAMPLE* » A* 

820 PRINT USING 830 * 'TENSILE STRENGTH (PSI) OF FIBER AT ’ » W1 r ’ GAGE LENG 
TH ISJ*»S 

830 IMAGE 34An4D.3Dn16AnD.3DE 

840 PRINT USING 850 i ’TOTAL NUMBER OF SPECIMENS TESTED ISJ’nNI 
850 IMAGE 42A . 2D 

860 PRINT USING 870 J ’TOTAL NUMBER OF BREAKS ISI’nN 
870 IMAGE 32A r 3D 
880 FOR 1=1 TO N 

890 PRINT USING 900 f ’D(’nIn’)n L(’nIn’)=’nD(I)nL<I) 

900 IMAGE 2An3Dn5An3Dn4An2D«2Dn5D.2D 

910 NEXT I 
920 GOTO 120 
930 ! 

940 ! EDIT ROUTINE ***** 

950 ! 

960 CLEAR 0 DISP * ADD ( A ) n CHANGE(C)» END(E)’J0 INPUT Z* 

970 IF Z#= ’ A ’ THEN 1000 
980 IF Z*=’C’ THEN 1110 
985 IF Z$= ’ E ’ THEN 120 
990 GOTO 960 
1000 ! ADD ROUTINE 

1010 DISP ’ENTER ADDITIONAL NUMBER OF SPECIMENS TESTED: ’J0 INPUT Ml 

1020 DISP ’ENTER ADDITIONAL NUMBER OF BREAKS « ’ » 0 INPUT M 

1030 N1=N1 +M1 

1040 FOR I=N+ 1 TO N+M 

1050 DISP ’D(’»If’)f L(’fl»’>=’ 

1070 INPUT D<I)rL(I) 

1080 NEXT I 

1090 N=N+M 

1100 GOTO 960 

1110 ! CHANGE ROUTINE 

1112 DISP ’CHANGE TENSILE STRENGTH (S>» MEAN CRITICAL LENGTH <W>» DIAMET 
ER AND LENGTH (L) ■ 

1113 INPUT D* 

1114 IF D*= ’S’ THEN 1117 

1115 IF D$= ’ L ’ THEN 1120 

1116 IF D$= ' W ’ THEN 1181 

1117 DISP ’TENSILE STRENGTHENS 

1118 DISP ’CHANGE TO TENSILE STRENGTH®’ 0 INPUT S 

1119 GOTO 1180 
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1120 DISP * ENTER ITEM * TO CHANGE' 

1130 INPUT I 

1140 DISP 'OLD VALUES:' 

1150 DISP "D('»I»")r L( ' f I» ' )=' »D(I) »L(I) 

1160 DISP 'NEW VALUES:' 

1170 DISP 'D('»I»')f L('flf*>='»0 INPUT D<I>»L(I) 

1180 GOTO 960 

1181 DISP 'MEAN CRITICAL LENGTH='»W1 

1182 DISP 'CHANGE TO MEAN CRITICAL LENGTH= ' 0 INPUT U1 

1183 GOTO 960 
1190 ! 

1200 ‘ *** UEIBULL DISTRIBUTION PARAMETERS *** 

1210 ! 

1220 SI r S2 r S3 * C » F2 = 0 

1230 ! INITIALIZE ALPHA TO BE 10 

1240 A= 1 0 

1250 ! INITIALIZE ITERATION STEP TO BE 1 
1260 X = 1 

1270 ! ASPECT RATIO DETERMINATION 
1280 FOR 1=1 TO N 
1290 R ( I ) =L < I ) /D ( I ) 

1295 NEXT I 

1299 CLEAR 

1300 DISP 'STORE DATA FOR HISTOGRAM /UEIBULL PLOT ? (Y OR N>'»0 INPUT Q* 
1310 IF Q*='Y' THEN 1330 

1320 IF Q$ = ' N ' THEN 1460 
1325 GOTO 1300 
1330 ! 

1340 ! STORE ASPECT RATIO FOR H I STO/ WE I BULL /NORMAL PLOT *** 

1350 ! 

1360 CLEAR 0 DISP 'ENTER FILE NAME? (END WITH H>'f0 INPUT B$ 

1370 DISP 'NEW OR EXISTING FILE (N/E)?' 

1380 INPUT P$0 IF UPC*(P$C1»1D)='N' THEN 1390 ELSE GOTO 1400 
1390 CREATE UPC*(B$>>20 
1400 ASSIGN* 1 TO UPC$(B*> 

1420 PRINT* 1 ? NfRO 
1440 ASSIGN* 1 TO * 

1450 DISP 'DATA ON 'JUPC*(B$)J' STORED ON TAPE.' 

1460 DISP 'DO YOU WANT A PRINT OUT? (Y/N)?' 0 INPUT C$ 

1461 IF C$ = * Y ' THEN 1463 

1462 IF C*= ' N ' THEN 1468 

1463 PRINT USING 1464 » *N'»'R' 

1464 IMAGE Ar8X»A 

1465 FOR 1=1 TO N 0 PRINT USING 1466 f I»R(I) 0 NEXT I 

1466 IMAGE 3D»5X»3D.3D 

1468 BEEP 

1469 CLEAR 

1470 SI » S2 » S3=0 
1480 FOR 1=1 TO N 

1490 S1=S1+R(I> ”A*LOG ( R ( I ) ) 

1500 S2=S2+R ( I ) ”A 
1510 S3=S3 4-LOG ( R ( I ) ) 

1520 NEXT I 

1525 DISP ' A= ' » A 

1526 DISP ' F2= ' ? F2 
1530 ! 

1540 F=S1/S2-1/A-S3/N 
1545 DISP " F= ' » F 
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1550 IF ABS<FX. 000001 THEN 1680 

1560 IF C>0 THEN 1600 

1570 F2=F 

1580 C=1 

1590 GOTO 1650 

1600 IF F2>0 THEN 1630 

1610 IF F<0 THEN 1650 

1620 GOTO 1640 

1630 IF F>0 THEN 1650 

1640 X=-X/10 

1650 A=A-X 

1660 F2=F 

1670 GOTO 1470 

1680 ! 

1690 B=(1/N*S2>~<1/A) 

1700 P = 1 - 1 / A 
1710 GOSUB 1960 
1900 G2=G 
1910 T=S/2/B*G2 

1920 P= 1 -2/A 

1921 DISP 

1930 GOSUB 1960 

1940 V=S"2/4/B~2*<G-G2”2> 

1941 U=SQR(V> 

1945 GOTO 2310 
1950 ! 

1960 ! GAMMA FUNCTION SUBROUTINE *** 

1970 ! 

1980 IF P - 5 7 < = 0 THEN 2020 

2000 G=1 ♦ E75 

2001 DISP * P>57 1 OVERFLOW t G SET TO l.E75'f 

2002 WAIT 3000 

2003 BEEP 
2010 RETURN 

2020 P1=P 
2030 E= « 00000 1 
2040 C=0 
2050 G= 1 

2060 IF Pl-2<=0 THEN 2110 

2070 IF Pl-2<=0 THEN 2220 

2080 P1=P1-1 

2090 G=G*P1 

2100 GOTO 2070 

2110 IF Pl-XO THEN 2140 

2120 IF Pl-1=0 THEN 2270 

2130 IF Pl-1>0 THEN 2220 

2140 IF PI -E>0 THEN 2190 

2150 Z=SGN(Pl)*INT(ABS(Pl)-.5)-Pl 

2160 IF ABS ( Z > -E<=0 THEN 2280 

2170 IF 1-Z-E<=0 THEN 2280 

2180 IF Pl-1>0 THEN 2220 

2190 G=G/P1 

2200 P1=P1+1 
2210 GOTO 2180 
2220 Z=P1-1 

2230 W=Z*<-.5684729+Z*< .2548205+Z*- .0514993 > ) 
2240 W=Z*< .985854+Z*(-.8764218+Z*< .8328212+U) ) ) 
2250 Gl=l+Z* ( - ♦ 5771017+W ) 



2260 G=G*G1 
2270 RETURN 

2280 DISP 'P IS WITHIN l.E-6 OF BEING A -VE INTEGER' 

2281 WAIT 3000 

2282 BEEP 
2290 RETURN 
2300 ! 

2310 ! DATA OUTPUT *** 

2320 PRINT USING 2330 i 'SAMPLE* ' f A$ 

2330 IMAGE 7Af10A 

2340 PRINT USING 2350 f 'TOTAL NUMBER OF SPECIMENS*' fNI 
2350 IMAGE 26Af3D 

2360 PRINT USING 2370 f 'TOTAL NUMBER OF BREAKS*' fN 
2370 IMAGE 23Af4D 

2380 PRINT USING 2390 i * WEIBULL DISTRIBUTION ALPHA='fA 
2390 IMAGE 27Af3D.3D 

2400 PRINT USING 2410 f 'WEIBULL DISTRIBUTION BETA='fB 
2410 IMAGE 26A r 4D ♦ 3D 

2420 PRINT USING 2430 f 'TENSILE STRENGTH (PSI) OF FIBER AT'fWIf' GAGE L 
ENGTH* ' t S 

2430 IMAGE 35A f 4D . 3D f 13A f D . 3DE 

2440 PRINT USING 2450 f 'INTERFACIAL SHEAR STRENGTH='fT 
2450 IMAGE 27AfD*3DE 

2460 PRINT USING 2470 i 'VARIANCE OF THE SHEAR STRENGTH^ fV 
2470 IMAGE 31AfD,3DE 
2480 ! 

2490 END 
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APPENDIX I 


WEI BULL DISTRIBUTION OF CRITICAL LENGTH 
A HISTOGRAM PLOT 


The purpose of this program is to plot the Weibull Distribu- 
tion curve and the population density curve of the critical aspect 
ratio of the fibers. 
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10 COM Z9 

20 ! ****#DATA PLOT PROGRAM***********VERSION 6 ♦ 09********* 

30 DIM P<2) ,L(2) ,N1<2) »S2<2) ,X (2,350) ,T*C803,X*C803,Y*C803,S*C803,L1<2) , 
Y1 < 2 , 350 ) fR<350) 

40 DIM T1*C803 

50 ! *****INITIALIZATION*********************************** 

60 G*= * AG * e N*=*NY * 8 K*='KTC* 8 L1<0)=0 8 S*=** 8 N=0 
70 N1=0 8 Z9= 1 8 DEG 
80 ON KEY# 1 , * ENTER ■ GOTO 180 
90 ON KEY# 2, ‘SETUP* GOTO 530 
100 ON KEY# 3, 'HEADING* GOTO 830 
110 ON KEY* 4 » * PLOT* GOTO 2130 
120 ON KEY* 5 , * STORE * GOTO 390 
130 ON KEY* 6t ’EDIT* GOTO 1460 
140 ON KEY# 7 » * LEGEND* GOTO 1390 
150 ON KEY* 8 » * RANGE* GOTO 900 
160 CLEAR 8 KEY LABEL 
170 GOTO 170 
180 ! 

190 ! *****DATA ENTRY ##*****#**#****###*#*#***#*#*##**#*** 

200 CLEAR 8 DISP *Enter Data from Keyboard(K) or Tape < T ) or from Calc u la 
tion(C) * 

210 INPUT Q$8 Q2=P0S(K*»UPC*(Q$C1»13) ) 

220 IF Q2=0 THEN 200 

230 ON 02 GOTO 240,750,741 

240 FOR I =N + 1 TO 6 

250 CLEAR 8 DISP 'Enter Data for Set *JI»* (Y/N)* 

260 INPUT Q*8 Q2=P0S<N*,UPC*<Q*C1,13) ) 

270 IF 02=0 THEN 250 
280 ON 02 GOTO 370»290 

290 DISP 'Enter Data Pai rs ( ENDLINE ends)* 

300 FOR J=1 TO 100 

310 DISP a X< a m a f a »J»*)v Y(*»If*f*»J»*) = * » 8 INPUT Q$,R* 

320 IF Q$= * * THEN 350 

330 X(I» J)=VAL(Q$) 8 Y1 ( I * J > =VAL ( R* ) 

340 NEXT J 

350 N 1 ( I ) = J-l 8 N1=N1+N1 ( I ) 

360 NEXT I 
370 N=I -1 

380 ON Z9 GOTO 530»160 

390 ! *****STORE DATA ON TAPE*** 

400 Z9=2 

410 CLEAR 8 DISP 'Store on Tape(Y/N) * 

420 GOSUB 3190 8 ON 09 GOTO 160,430 
430 DISP 'FILE NAME * J 8 INPUT F* 

440 DISP * New or Existing File (N/E)* 

450 INPUT Q$8 IF UPC* ( Q*C 1 , 1 3 > = * N * THEN 460 ELSE GOTO 470 
460 CREATE UPC*<F*)»50 
470 ASSIGN* 1 TO UPC*(F*) 

480 PRINT# 1 i N f N1 ( ) , XO , YO ,X1 , Y1 » X ( , > , Y1 ( , ) , LI ( ) , P ( ) , L ( > , S ( ) , S* , T* » Y* , X 

*fF1fF2fQ1fX2fX3fY2fY3 

490 PRINT* 1 * Tl* 

500 ASSIGN# 1 TO * 

505 CLEAR 

510 DISP * DATA SET *fUPC*(F*)f* STORED ON TAPE* 8 KEY LABEL 
520 Z9=2 8 GOTO 170 

530 ! *****DETERMINE COLOR, LINE AND SYMBOL************** 

540 CLEAR 
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550 FOR 1=1 TO N 
560 DISP 

570 DISP 'ENTER PEN * <l-4>»' 0 DISP 'LINETYPE # (1-8) AND PL 

OTTING SYMBOL * <l-6>' 

580 DISP "FOR DATA SET'JI 
590 IF 1=1 THEN 600 ELSE 700 

600 PRINT 'PEN COLOR KEY ■ 0 PRINT '1= BLACK' 0 PRINT *2 RED' 0 PRINT '3 
GREEN' 0 PRINT '4 BLUE' 

610 PRINT 0 PRINT 'LINETYPE KEY' 0 PRINT '1= SOLID' 

620 PRINT '2 BLANK' 0 PRINT '3 DOT' 0 PRINT '4 S.DASH' 

630 PRINT *5 L » DASH * 0 PRINT '6 DASH/DOT' 

640 PRINT *7 L ♦ DASH/S . DASH ' 0 PRINT '8 DASH/DOT/DOT' 

650 PRINT 0 PRINT 

660 PRINT 'PLOTTING SYMBOL KEY' 

670 PRINT '1= DOT' 0 PRINT '2 CIRCLE' 0 PRINT '3 SQUARE' 

680 PRINT '4 TRIANGLE* 0 PRINT *5 PLUS SIGN' 0 PRINT *6 DIAMOND' 

690 PRINT 0 PRINT 0 PRINT 0 PRINT 
700 INPUT P(I)»L(I)fS(I) 

710 IF P(I)<1 OR P ( I > >4 THEN 700 

720 IF L< I XI OR L ( I ) >8 THEN 700 

730 IF 8(1X1 OR S ( I ) >6 THEN 700 

740 NEXT I 0 ON Z9 GOTO 830»160 

741 ! CALCULATION FOR UEIBULL DISTRIBUTION 

742 DISP 'ENTER FILE NAME WITH R() MATRIX FOR UEIBULL DISTRIBUTION* 0 IN 
PUT F*0 ASSIGN* 1 TO F* 

743 READ* 1 t N5»R< ) 

744 ASSIGN* 1 TO * 

745 N= 1 0 N1 < 1 ) =N5 

746 DISP 'ENTER VALUES' 0 INPUT A » B 

747 CLEAR 0 DISP 'CUMULATIVE PLOT (C) OR DENSITY POPULATION PLOT <D>?' 0 
INPUT E* 

748 IF E$= ' C ' THEN 750 

749 IF E* = ' D ' THEN 752 ELSE 747 

750 FOR 1 = 1 TO N5 0 X < 1 » I ) = ( R ( I ) /B ) "A 0 Y1 ( 1 » I > = 1-EXP ( -X < 1 » I > ) 0 X(1»I) = 
R ( I ) 0 NEXT I 

751 Z9=2 0 GOTO 160 

752 FOR 1=1 TO N5 0 X ( 1 > I >= ( R ( I ) /B> "A 0 Y1 ( 1 , I ) =A/B* ( R ( I ) /B ) " < A-l ) *EXP ( - 
X < 1 f I ) ) 0 X < 1 » I ) =R < I ) 

753 NEXT I 

754 Z9=2 0 GOTO 160 

759 ! DATA ENTRY*************** 

760 DISP 'ENTER FILE NAME ' i 0 INPUT F* 

770 ASSIGN* 1 TO UPC*<F*> 

780 READ* 1 i N»Nl()»XO»YOfXlfYl»X(»)fYl(»)»Ll()»P()»L()»S()f S$ » T$ » Y* » X* 
fF1fF2fQ1» X2fX3fY2fY3 
790 READ* 1 i Tl* 

800 ASSIGN* 1 TO * 

810 Z9=2 
820 GOTO 160 

830 ! ##**#ENTER HEADINGS****************************** 

840 CLEAR 0 DISP 'ENTER TITLE' *0 INPUT T* 

850 DISP 'ENTER SUBTITLE* *f0 INPUT T1S 
860 DISP 'ENTER X axis TITLE' »0 INPUT X* 

870 DISP 'ENTER Y axis TITLE' f0 INPUT Y* 

880 P<0)=1 

890 ON 19 GOTO 900*160 

900 ! *****DEFINE GRAPH******************************** 

910 XO t Y0 = 10"99 0 XI r Yl = -10"99 
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1010 FOR 1=1 TO N 

1020 FOR J=1 TO N1(I> 

1030 XO=MIN(XO*X(I* J> ) B X1=MAX < XI *X < I » J) > 

1040 Y0=MIN(Y0*Y1(I*J>) B Y1=MAX (YliYKIr J)) 

1050 NEXT J 
1060 NEXT I 

1070 DISP ’X DATA RANGE IS ’JXO*’ TO ’JX1 

1080 DISP * Y DATA RANGE IS ’*Y0*’ TO ’*Y1 

1100 DISP * DO YOU WISH TO SET GRAPH LIMITS ’*8 GOSUB 3170 

1110 ON Q9+1 GOTO 1 100* 1 160 * 1 120 

1120 DISP ’ ENTER Xmin*Xm3X* Yain»Y*ax ’» 

1140 INPUT X0*X1*Y0*Y1 

1150 IF X1<X0 OR Y1<Y0 THEN 1120 

1160 DISP 'DO YOU WISH AXES OR GRID <A or G)’* 

1170 INPUT 0*8 Q1=P0S < G* » 0* ) 

1175 DISP B DISP ’NUMERICAL LABELS* ’ 

1177 DISP ’on horizontal axis...’ 

1180 ON Ql+1 GOTO 1 160* 1 190 * 1240 

1190 DISP ’Units per label* units per tic mark’? 

1200 INPUT X3 * X2B X3=INT < X3/X2+ ♦ 51 ) B IF X3<0 THEN 1190 
1205 DISP B DISP ’on vertical axis..*’ 

1210 DISP ’Units per label* units per tic mark.’? 

1220 INPUT Y3 » Y2S Y3=INT ( Y3/Y2 ) 8 IF Y3<0 THEN 1210 
1230 GOTO 1290 

1240 DISP ’units per GRID LINE*units per tic mark’* 

1250 INPUT X3 * X26 X3=INT ( X3/X2+ . 51 ) 8 IF X3<0 THEN 1240 
1255 DISP B DISP ’on vertical axis...’ 

1260 DISP ’units per GRID LINE*units per GRID LINE'* 

1270 INPUT Y3*Y2B Y3=INT < Y3/Y2+ .51 > 8 IF Y3<0 THEN 1260 

1280 ON ERROR GOSUB 3220 

1290 IF Y2<1 THEN F2=2 ELSE F2=0 

1300 IF Y2C.1 THEN F2=4 

1310 IF X2<1 THEN Fl=2 ELSE F1=0 

1320 IF X2C.1 THEN Fl=4 

1330 DISP ’INTEGER LABELS ON X-AXIS’ » 

1340 GOSUB 3170 

1350 ON Q9 GOTO 1370*1360 

1360 F1=0 

1370 OFF ERROR 

1380 ON Z9 GOTO 1390*160 

1390 ! *******ENTER LEGENDS****** 

1400 S*=’ ’ 8 CLEAR 
1410 FOR 1=1 TO N 

1420 DISP ’Enter Lesend for Group '*1 
1430 INPUT I*B S*=S*XI* 8 LI ( I ) =LEN ( S* ) 

1440 NEXT I 

1450 ON Z9 GOTO 390*160 
1460 ! EDIT ROUTINE 

1470 CLEAR 8 DISP ’Enter data set number’ 

1480 INPUT I 

1490 DISP ’There are’JNKI)?’ data pairs <X*Y) in this set’ 

1500 DISP 8 DISP ’Enter startinS % endind item numbers for list (ENDLINE 

ends) ’ 

1510 INPUT D* * E* 

1520 IF D*=” THEN 2122 

1530 FOR J=VAL ( D* ) TO MIN ( UAL (E* ) »N1 < I ) ) 

1540 DISP J*’=’*X(I»J)*’*’*Y1(I*J) 

1550 NEXT J 
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1560 DISP ' Add( A) » Delete(D) »Ch3n«le(C> * Insert < I ) rEnd<E> ■ 

1570 INPUT Yl* 

1580 IF UPC$( YliCl il])='A' THEN 1640 
1590 IF UPC$<Yl*Clf 1])='D' THEN 1700 
1600 IF UPC$<Yl*Clrl3>=*C* THEN 2050 
1610 IF UPC*<Y1*C1» 13>=*I* THEN 1850 
1620 IF UPC*<Yl*Cl»n> = *E* THEN 1500 
1630 GOTO 1500 
1640 ! ADD ROUTINE 

1650 DISP , X( , }I! , » , ;Nl(I)+l» , )fY(SIJ , *'fNlU)H»') = * 

1660 N1 < I )=N1 ( I )+l 

1670 INPUT X(I>N1(I))»Y1(I»N1(I1) 

1680 GOTO 1560 
1690 ! 

1700 ! DELETE ROUTINE 

1710 DISP * Enter item # to delete* 

1720 INPUT J 

1730 DISP *X<* JIJ*r*f JJ *>rY(*JIJ*»*JJJ*> = *JX<I»J>J* * v Yl < I f J) 

1740 DISP * Enter D to delete* 

1750 ! 

1760 INPUT Yl* 

1770 IF UPC*< Y1*C1 »13>**D* THEN 1560 
1780 N1 ( I ) =N1 < I ) -1 
1790 FOR J1=J TO N1 ( I ) 

1800 X(I» J1)=X(I» Jl+1) 

1810 ! 

1820 Y1<I» Jl)=Yl(If Jl+l) 

1830 NEXT J1 

1840 DISP •Item* » J» * deleted* 0 GOTO 1560 
1850 ! INSERT ROUTINE 

1860 DISP 'Insert item before item ** 

1870 ! 

1880 INPUT J 

1890 DISP *X(*Ji;*f *J J» *)fY(*>IJ*f *J JJ’> = * 

1900 INPUT G8 f G9 

1910 N1 ( I ) =N1 ( I > +1 

1920 FOR J1=N1 ( I ) TO J+l STEP -1 

1930 ! 

1940 X(I»J1)>X(I»J1-1> 

1950 Y1(I»J1)-Y1(I»J1-1> 

1960 NEXT J1 
1970 X ( I » J ) =G8 
1980 Yl ( I » J ) =G9 
1990 ! 

2000 GOTO 1560 
2010 ! 

2020 ! 

2030 ! 

2040 ! 

2050 ! CHANGE ROUTINE 

2060 DISP * Enter Item # to Change* 

2070 INPUT J 

2080 DISP * Old Values! ' 0 DISP 
2090 DISP X(I»J)»*»*»Yl(IrJ) 

2100 DISP 0 DISP 'Enter New Values* 

2110 INPUT X<I»J)»Y1(I»J)0 DISP 0 DISP 
2120 GOTO 1560 

2122 CLEAR 0 PRINT USING *9/* » 0 PRINT 'Press RANGE User Kes before 
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PLOT ■ 0 KEY LABEL 0 GOTO 170 

2130 ) *****GRAPH****************************************** 

2140 DISP * PLOT ON ♦ * 

2150 DISP ■ 1. HP SCREEN* 

2160 DISP * 2. PLOTTER* 0 DISP * 3« POSTER* 

2170 INPUT Z 

2180 IF Z=1 THEN 2270 

2190 IF Z=2 THEN 2230 

2200 IF Z=3 THEN 2230 

2210 GOTO 2140 

2220 Z9=2 

2230 PLOTTER IS 705 
2240 IF Z=3 THEN 2280 
2250 LIMIT 0 f 250 f 15 f 195 
2260 GOTO 2280 
2270 PLOTTER IS 1 
2280 LOCATE 26f120f20f80 
2290 SCALE XOfXIfYOfYI 

2300 FXD FI f F2 0 CSIZE 4f»6f0 0 LINETYPE 1 0 GCLEAR 0 PEN 1 
2310 ON Q 1 GOTO 2320f2340 
2320 LAXES -X2 f Y2 f XO f YO » X3 f Y3 
2330 GOTO 2350 

2340 LGRID -X2 f Y2 f XO f YO f X3 f Y3 0 FRAME 

2350 PEN P(0) 0 MOVE X0+ ( Xl-XO > /2 f Y 1+ < Yl- YO ) *3/20 

2360 LORG 4 0 CSIZE 2*Z+2f,7fO 

2370 LABEL T* 

2380 CSIZE 2*Z + 1 f , 7 
2390 LABEL Tl$ 

2400 SETGU 
2410 MOVE 70 f 3 
2420 LORG 4 
2430 LABEL X$ 

2440 MOVE 8 f 45 

2450 LDIR 90 0 LABEL Y* 

2460 LDIR 0 
2470 SETUU 

2650 ! *****PLOT DATA************************************ 

2660 FOR 1=1 TO N 
2670 PEN P<I> 

2680 MOVE X(IfI) fYKIfI) 

2690 GOSUB 2890 

2700 MOVE X(IfI)fYKM) 

2710 FOR J=2 TO N1(I) 

2720 LINETYPE L<I> 

2730 DRAW X ( I f J ) f Y l < I f J ) 

2740 GOSUB 2890 

2750 MOVE X(IfJ)fYKIfJ) 

2760 NEXT J 
2770 GOSUB 2800 
2780 NEXT I 0 PEN 0 

2782 IF Z=1 THEN CLEAR 0 DISP USING *7/* i 0 DISP 'Press KEY LABEL when 
ready* 0 WAIT 2000 

2784 IF Z=1 THEN CLEAR 0 GRAPH 0 GOTO 170 
2790' GOTO 160 

2800 ! *****DRAW KEY************************************** 

2810 SETGU 0 LINETYPE L<I> 

2820 MOVE 1 10 f 96-3* ( I -1 ) 

2830 IDRAW 5f0 
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2840 GDSUB 2890 0 SETGU 
2850 MOVE 1 17 f 96-3# ( I-i > 

2860 LORG 2 8 CSIZE 3f.6 0 LOIR 0 
2870 LABEL S*CL1 ( 1-1 ) +1 f LI ( I ) 3 
2880 SETUU 8 RETURN 

2890 ! *****PLOT SYMBOLS********************************** 
2900 SETGU 0 LINETYPE 1 

2910 ON 3<I> GOTO 2920, 2940 f 2990 f 3040 f 3080 f 3130 
2920 ! ##*##DOT#*#**#*#***#**##* 

2930 RPLOT OfO 8 SETUU 0 RETURN 
2940 ! *****CIRCLE************** 

2950 IMOVE *15*0 

2960 FOR A=0 TO 360 STEP 30 

2970 PDIR A 0 RPLOT *8f0 

2980 NEXT A 0 SETUU 0 RETURN 

2990 ! **** ^SQUARE ******** ****** 

3000 IMOVE -1 f -1 

3010 IDRAW 2 f 0 0 IDRAW 0f2 

3020 IDRAW -2f0 0 IDRAW Of-2 

3030 SETUU 0 RETURN 

3040 ! **#**TRIANGLE************ 

3050 IMOVE OfI 

3060 IDRAW 1 f -2 0 IDRAW -2,0 
3070 IDRAW 1 f 2 8 SETUU 0 RETURN 
3080 ! *****PLUS SIGN*********** 

3090 IDRAW OfI 0 IDRAW Of-2 

3100 IMOVE OfI 

3110 IDRAW 1 f 0 0 IDRAW -2f0 

3120 SETUU 0 RETURN 

3130 ! * #***D I AMOND* ************ 

3140 IMOVE OfI 0 IDRAW 1f-1 
3150 IDRAW -1 f -1 0 IDRAW -1f1 
3160 IDRAW 1 f 1 0 SETUU 0 RETURN 

3170 ! *****YES OR NO************************************** 

3180 DISP ' (Y/N> * 

3190 INPUT Q*0 Q9=P0S(N*fQ*) 

3200 IF Q9=0 THEN 3180 
3210 RETURN 

3220 ». *****CHECK FOR DIVISION BY************ZERO******* 
3230 IF ERRN=8 THEN FI f F2 = 0 
3240 RETURN 



APPENDIX J. INSTRUMENTED IMPACT TESTER 


The ETI-Dynatup drop dart instrumented impact tester monitors 
and records the entire impact event, from initial impact and 
acceleration from rest to plastic bending to fracture initiation 
and propagation to failure. Also, the drop dart system permits 
dynamic tests to be made at impact energies and velocities not 
available with other testing machines. 

The drop dart impact test consists of a drop tower and instru- 
mentation package - a microprocessor. The drop tower is a gravity 
driven impact machine equipped with remote controls for release 
of the hammer and tup assembly, and a motorized lift mechanism 
for easy return of the hammer to a predetermined drop position. 

An automatic rebound brake is furnished as standard. 

The hammer and tup assembly is designed so that the weight of 
the hammer can be varied and the tups can be interchanged. The 
base plate of the drop tower permits interchangeability of anvil 
supports to accommodate different specimen size and geometries. 

A typical output is shown in Figure J-l. Pj and Pp are the loads 
correspond to the initial damage of the composite and the maximum 
load required for punch through. The corresponding energy ab- 
sorbed at each stage is Ej and Ep. The total energy absorbed 
during the whole event is E T . In this report, Pp and E T are used 
to characterize the impact toughness of the laminates. 


J-l 



Load— N (lb) 


1792 

(400) 


1344 Pi 
(300) f 1 


896 

( 200 ) 


448 

(100) 
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